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ABSTRACT 

Context. Near- to mid-infrared observations of molecular emission from protoplanetary disks show that the inner regions are rich in 
small organic volatiles (e.g., C^H 2 and HCN). Trends in the data suggest that disks around cooler stars (T eff ~ 3000 K) are potentially 
(i) more carbon-rich; and (ii) more molecule-rich than their hotter counterparts (T e fj > 4000 K). 

Aims. To explore the chemical composition of the planet-forming region (< 10AU) of protoplanetary disks around stars over a range 
of spectral types (from M dwarf to Herbig Ae) and compare with the observed trends. 

Methods. Self-consistent models of the physical structure of a protoplanetary disk around stars of different spectral types are coupled 
with a comprehensive gas-grain chemical network to map the molecular abundances in the planet-forming zone. The effects of 
(i) N, self shielding; (ii) X-ray-induced chemistry; and (iii) initial abundances, are investigated. The chemical composition in the 
‘observable’ atmosphere is compared with that in the disk midplane where the bulk of the planet-building reservoir resides. 

Results. M dwarf disk atmospheres are relatively more molecule rich than those for T Tauri or Herbig Ae disks. The weak far-UV flux 
helps retain this complexity which is enhanced by X-ray-induced ion-molecule chemistry. N 2 self shielding has only a small effect in 
the disk molecular layer and does not explain the higher C 2 H,/HCN ratios observed towards cooler stars. The models underproduce 
the OH/HtO column density ratios constrained in Herbig Ae disks, despite reproducing (within an order of magnitude) the absolute 
value for OH: the inclusion of self shielding for H.,0 photodissociation only increases this discrepancy. One possible explanation is 
the adopted disk structure. Alternatively, the ‘hot’ H,0 (T > 300 K) chemistry may be more complex than assumed. The results for 
the atmosphere are independent of the assumed initial abundances; however, the composition of the disk midplane is sensitive to the 
initial main elemental reservoirs. The models show that the gas in the inner disk is generally more carbon rich than the midplane ices. 
This effect is most significant for disks around cooler stars. Furthermore, the atmospheric C/O ratio appears larger than it actually is 
when calculated using observable tracers only. This is because gas-phase CL is predicted to be a significant reservoir of atmospheric 
oxygen. 

Conclusions. The models suggest that the gas in the inner regions of disks around cooler stars is more carbon rich; however, calcula¬ 
tions of the molecular emission are necessary to definitively confirm whether the chemical trends reproduce the observed trends. 
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1. Introduction 

Protoplanetary disks provide the ingredients - dust, gas, and ice 
- for planets and planetesimals such as comets (for a review, see, 
.g., | Williams & Cieza]|201 1| >. In disks around low-mass stars 
(< 2 M 0 ), planetary systems are thought to form relatively close 
to the parent star (< 10 AU); hence, the chemical composition 
of the inner disk region sets the initial conditions and avail¬ 
able elemental components for planetary systems. The molecu¬ 
lar material within ~ 10 AU is generally dense (> 10 8 cm~ 3 ) and 
can reach high temperatures (> 100 K) allowing (ro)vibrational 
excitation of molecules which emit radiation at near- to mid- 
infrared (IR) wavelengths. The physical conditions within this 
region drive the chemistry towards the formation of small, sim¬ 
ple, stable molecules most of which, fortunately, also have strong 
(ro)vibrational transitions. 

Near- to mid-IR spectroscopy of nearby primordial proto¬ 
planetary disks has demonstrated that the inner planet-forming 
regions are rich in organic volatiles. The Spitzer Space Telescope 


allowed the first detection of simple organic molecules in proto¬ 
planetary disks at IR wavelengths. Lahuis et al.|( j2006| ) reported 
absorption bands of C 2 H 2 , HCN, and C0 2 in the spectrum of the 
low-mass young stellar object, IRS 46, attributed to absorption 
from hot molecules in a disk within a few AU of the embedded 
star. The following year, |Gibb et akj ( |2QQ7[ > detected absorption 
bands from CO, C 2 H-,, and HCN in the disk around a member 
of the binary system, GV Tau, using NIRSPEC on Keck. Ongo¬ 
ing efforts have detected additional molecules in either emission 
or absorption, including, OH, H-,0, and CH 4 , in several nearby 
disks using both Spitzer and ground-based facilities (see, e.g., 


Carr & Najita 

2008; Salyk etal.J2008 

|Pascucci et al.|2009; Pon- 

toppidan et a 

. 2010;lCarr & Najita 

2011 Fedele et al. 2011; 

Salyk etal. 2011 

Mandell et al. 2012 

Bast et al. 20131 Najita et 

al. 2013; G bb & 

Torne 2013. Pascueci et al.12013). Because of 


dust opacity, such observations probe the composition of the disk 
atmosphere only. It remains unclear whether the atmospheric 
composition is representative of that of the disk midplane within 
which planetesimals sweep up the bulk of their material. 
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Several interesting trends have been noticed in the IR ob¬ 
servations. Pascucci et al. (2008 2009) presented results from a 
low-resolution (R as 64 - 128) Spitzer/IRS survey of more than 60 
sources ranging from brown dwarfs ( T e ff ~ 3000 K) to Sun-like 
stars CT eff ~ 5000 K). The observations demonstrated an under¬ 
abundance of HCN relative to C 2 H 2 in disks around M dwarfs 
compared with those around T Tauri stars. The authors postulate 
this could be due to the difference in FUV (far-ultraviolet) lu¬ 
minosity: M dwarf stars have insufficient FUV photons to disso¬ 
ciate N 2 , the main nitrogen reservoir, thereby trapping elemen¬ 
tal nitrogen which would otherwise be available to form other 
nitrogen-containing species, e.g., HCN. The authors have since 
published additional observations at higher spectral resolution 
(.R ~ 600) and find the same result: the ratio of C 2 H,/HCN line 
flux and relative column density decreases with increasing spec- 


10 - 36 pm wavelength range, with a source list also covering 
Herbig Ae/Be stars (T e s > 10,000 K). They detected strong H 2 0 
mid-IR line emission from 22 T Tauri stars in their sample, with 
a detection rate on the order of 2/3; however, to their surprise, no 
disks in their sample of 25 Herbig stars exhibited water (nor OH) 
line emission. At near-IR wavelengths, Fedele et al. ( 2011) con¬ 
ducted a high-resolution spectroscopic survey (L-band) of Her¬ 
big Ae/Be disks. OH emission was detected in only four objects, 
mainly flared disks, and, similar to that found by |Pontoppidan 


et al. (2010 1 , H 2 0 was not detected. Both sets of authors sug- 


fturry of astrochemical models mainly focussed on T Tauri disks 


(e.g, Agundez et al. 2008 

Willacy & Woods 2009; 

Woods & 

Willacy 2009; Walsh et al. 

2010,(Najita et al. 2011). T 

tese mod- 


[ak| ( |2008j ) adapted a model used for photon-dominated regions 
(PDRs) but neglected heating due to UV excess emission and X- 
rays and also assumed that the dust and gas temperatures were 
equal. |Willacy & Woods| ( |2009| > and |Woods & Willacy| ( |2009| > 
calculated the gas temperature by solving the equation of ther¬ 
mal balance; however, the dust temperature and density were as¬ 
sumed and heating by UV excess emission from the star was 
neglected ( D’Alessio et al.|2006| ). Najita et al. ( j2011 j > adopted a 


similar approach albeit using a reduced chemical network. Walsh 


|et al.| ( |20101 > used a self-consistent model for the protoplane¬ 
tary disk structure including the effects of heating by UV ex 
cess emission and X-rays ((Nomura & Millar|j2005[ [Nomura et 


|aLl|2007) and a chemical network similar in complexity to the 


work of [Willacy & Woods ( 2009j >, excluding deuterium. Sev 


eral works have also concentrated solely on water production in 


tral type (Pascucci et al. 2013 i. They also find that the HCN/H 2 0 
line flux ratios are higher for M dwarf stars than for T Tauri stars 
leading the authors to conclude that the C/O ratio in the inner 
regions of M dwarf and brown dwarf disks is higher (~ 1) than 
that for disks around T Tauri stars (< 1). Within the subset of 
T Tauri disks, Najita et al. ( 2013] ) postulated that a second trend 
was present. They found a correlation between the HCN/H 2 0 
line flux ratio and disk mass. Their hypothesis is that planetesi- 
mal formation is more efficient in higher mass disks and is able 
to lock up a significant fraction of oxygen (in the form of wa¬ 
ter ice) thereby increasing the C/O ratio in the inner region of 
the disk. In summary, the C/O ratio in the inner regions of proto¬ 
planetary disks appears to increase with decreasing spectral type, 
and within a particular sub class of star-disk systems, to increase 
with increasing disk mass. 

Pontoppidan et al. ( 2010j > conducted a similar survey in the 


the inner regions of protoplanetary disks (Glassgold et al. 

2009 

Bethell & Bergin 2009; V 

eijerink et al. 20091 Woitke et al. 

2009 

Adamkovics et al. 2014-j 

Du & Bergin 2014); however, to date. 

there has been no detailec 

study of chemistry in the inner regions 


of M dwarf nor Herbig Ae stars, especially to address the trends 
seen in the C 0 H 2 /HCN line ratios. 

In this work, we compute the physical and chemical structure 
of the planet-forming region (< 10 AU) of a protoplanetary disk 
around stars of different spectral types: (i) an M dwarf star; (ii) a 
T Tauri star; and (ii) a Herbig Ae star. Our aim is to investigate 
whether the stellar radiation field plays a role in the observed 
trends in CiH,,/HCN and OH/H 2 0 ratios derived from the mid- 
IR observations. Given the proposed importance of N 2 photodis¬ 
sociation for the formation of HCN, we investigate the effect of 
N 2 self shielding using recently computed shielding functions 


gested that the stronger FUV flux from the Herbig Ae/Be stars 
dissociates molecules in the unshielded inner disk region thereby 
lowering the line flux in relation to the stellar luminosity. Obser¬ 
vations from Herschel support this hypothesis: most of the Her¬ 
big Ae/Be disks in the GASPS and DIGIT key programs exhibit 
strong OI emission 63 pm with a sub sample also showing OH 
emission (Mee us et al.|2012] [Fedele et al.|2013] ). Many of these 
sources also have strong OI emission at 6300 A (see e.g., Acke 
|et a1. 1 2 00? i. In contrast, only HD 163296 (which is a ‘flat’ or set¬ 
tled disk) has a robust (> 3 cr) water detection between 50 and 
220 pm ( |Fedele et al .|20T2] [20T3] >. 

One of the first chemical models which concentrated solely 
on the inner disk (< 10 AU) of a T Tauri star is presented in 
Mar kwick et al.| (2002); however, the authors neglected the in¬ 
fluence of the stellar and interstellar radiation fields on the disk 
physics and chemistry using the reasoning that viscous heating 
will dominate the inner disk structure. The initial detections of 
volatiles in the inner regions of protoplanetary disks ([Lahuis 
|et al.||2006( Gibb et ak||2007| Carr & Najita||2008 1 prompted 


(Li et al. 2013a[ > on the subsequent nitrogen chemistry. We also 
use these models to probe the connection between the ‘observ¬ 
able’ gas emission from the disk atmosphere (at near- to mid- 
IR wavelengths) with the chemical composition of the midplane 
within which forming planets and planetesimals sweep up the 
bulk of their elemental building blocks. 

The remainder of the paper is stmctured as described. In 
Sect. [2] we outline the methods for computing the disk physi¬ 
cal and chemical structure, in Sect.[3]we present our results, and 
in Sects. [4] and [5] we discuss the implications of this work and 
summarise the main results, respectively. 


2. Protoplanetary disk models 

2.1. Physical model 

The physical structure of each disk model is calculated using the 
methods outlined in [Nomura & Millar ( 2005) w ith the addition 
of X-ray heating as described in Nomura et al. ( 2007) . Because 
the methodology is covered in detail in a series of previous pub¬ 
lications (see, e.g., Nomura & Millar|2005|[Nomura et akj2007 ; 


els differed somewhat in their level of complexity. |Agundez~et| 
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Walsh et al.|2010 |2012 2014) , here, we highlight the important 
parameters only. 

For each disk model, we assume the disk is steady, ax- 
isymmetric, and in Keplerian rotation about the central star. We 
parametrise the kinematic viscosity via the dimensionless a pa¬ 
rameter which scales the maximum size of turbulent eddies by 
the disk scale height, //, and the sound speed of the gas, c s , i.e., 
v ~ aHc s . For protoplanetary disks, a ~ 0.01. We model the 
structure of a disk surrounding a star of three different spectral 
types: an M dwarf star, a T Tauri star, and a Herbig Ae star. 
We list the stellar mass, M*, stellar radius, R*, and effective 
temperature, R*, of each host star in Table [T] We also list the 
adopted mass accretion rate of each star-disk system, M , and 
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the gas mass surface density at 10 AU, Sioau- We assume the 
T Tauri and Herbig Ae systems have a mass accretion rate typ¬ 
ical for pre-main-sequence stars, ~ 10 s M 0 yr 1 . Because the 
accretion signatures from lower-mass stars are not as strong as 
for the higher-mass systems, we assume a lower accretion rate 
for the M dwarf system, ~ 10 9 M G yr -1 (see, e.g., Herczeg 
|& Hinenbrand|2009l ). We calculate the dust temperature assum¬ 
ing local radiative equilibrium between the absorption and ree¬ 
mission of radiation by dust grains. For the calculation of the 
FUV extinction by dust grains, we adopt a dust-grain size dis¬ 
tribution which replicates the extinction curve observed in dense 
clouds (Weingartner & Drain epOOl [ i. We calculate the gas tem¬ 
perature assuming detailed thermal balance between the heat¬ 
ing and cooling of the gas. We include heating via photoelectric 
emission from dust grains induced by FUV photons and heat¬ 
ing due to the X-ray ionisation of H atoms. The cooling mecha¬ 
nisms included are gas-grain collisions and line transitions. The 
radiation field of each star is simulated as a black body at the 
stellar effective temperature (as listed in Table [lj with UV ex¬ 
cess emission scaled to the relative mass accretion rates for the 
M dwarf and T Tauri disks. The stellar FUV (912 - 2100 A) 
radiation field at 1 AU is shown in Figure [T] for all three cen¬ 
tral stars. The UV excess emission has two components derived 
from a best-fit model of the observed TW Hya spectrum: a di¬ 
luted black-body spectrum to simulate bremsstrahlung emission 
(7br ~ 25,000 K) and Lyman-cr line emission. The Lyman-a 
line is modelled as a Gaussian with a width « 2 A and the peak 
flux is determined assu ming a line/continuum luminosity ratio 
of 10 3 (see Nomura & Millar 2005 and references therein). Be- 


and see Sect. 2.2.5[ ). We have also added a set of three-body re 
actions compiled for use in combustion chemistry model^](see. 


body processes become increasingly important in the inner disk 
region where the density and temperature are sufficiently high 
(> 10 10 cm -3 and > 1000 K). We have also included at least one 
collisional dissociation reaction (AB + M —> A + B + M) for all 
neutral species expected to be abundant in the inner disk. We in¬ 
clude a small chemical network involving vibrationally excited, 
or ‘hot’, H,. For each gas-phase neutral-neutral reaction involv¬ 
ing FL which also possesses an activation barrier, we include a 
duplicate reaction involving hot H 2 with a barrier reduced by the 
internal energy of the excited H , (a; 30,163 K, see, e.g., Bruderer 
|et al.|2012] and references therein). 

2.2.2. Gas-grain interactions 

We allow the freezeout (adsorption) of molecules on dust grains 
forming ice mantles and the desorption (sublimation) of ices via 
thermal desorption and photodesorption (|Tielens & Hagen 1982 


cause young stars also exhibit strong X-ray emission, we use a 
TW Hya-like X-ray spectrum (generated by fitting the observed 
XMM-Newton spectrum) for the M dwarf and T Tauri stars, with 
a total luminosity, L x ~ 10 30 erg s' 1 (Preibisch et al.]2005) and 
assume L x ~ 3 x 10 29 erg s 1 and T x ~ 1.0 keV for the X-ray 
spectrum of the Herbig Ae star (see, e.g., Zinnecker & Preibisch 
|1994||Hamaguchi et al.|2005] l. 

2.2. Chemical model 

The network we use to calculate the disk chemical evolution in¬ 
cludes gas-phase reactions, gas-grain interactions (freezeout and 
desorption), and grain-surface chemistry. 

2.2.1. Gas-phase network 

The basis for the gas-phase chemistry is the complete net¬ 
work from the recent release of the UM1ST Database for As- 
trochemistry (UDfA) termed ‘Rate12’ which is publicly avail¬ 
ably (McElroy et al. 2013 1 . Rate12 includes gas-phase two- 
body reactions, photodissociation and photoionisation, direct 
cosmic-ray ionisation, and cosmic-ray-induced photodissocia¬ 
tion and ionisation. In this work, the photodissociation and pho¬ 
toionisation rates are calculated by integrating over the specific 
reaction cross section for each species and the calculated FUV 
spectrum at each p oint in the disk (|Walsh et al.|2012 1 using the 
cross sections from |van Dishoeck et al. ( 2006]) . Similar to pre¬ 
vious work, we have supplemented this gas-phase network with 
direct X-ray ionisation reactions and X-ray-induced ionisation 
and dissociation processes (as described in [Walsh et al.|[2012 


Hasegawa et al.|1992| Walsh et al.|2010 2012| >. We adopt the set 

of molecular binding energies compiled for use in conjunction 
with Rate12 ( McElroy et al.j20T3) >. We have updated the bind¬ 
ing energies in light of recent measurements for HCN ( [Noble et] 
[ak||2012j ). To simplify the calculation of the gas-grain interac¬ 
tion rates, we assume compact spherical grains with a radius of 
0.1 /rm and a fixed density of ~ 10 12 relative to the gas num¬ 
ber density. Each grain thus has ~ 10 6 surface binding sites. We 
include photodesorption by both external photons and photons 
generated internally via the interaction of cosmic rays with H 2 
molecules. We use experimentally determined photodesorption 
yields where available (see, e.g., Oberg et al. 2009abe). For 
all other species we use a yield of 10 1 molecules photon -1 . 
In the calculation of the freezeout rates, we assume a sticking 
coefficient, S ~ 1, for all species except H, for which we use a 
temperature-dependent expression which takes into account both 
physisorption and chemisorption and describes the decreased 
sticking probability at higher temperatures (Sha et al. 2005 [ Cup 


pen et al. 2010b I. We assume the rate of H 0 formation equates to 
half the rate of arrival of H atoms on dust grain surfaces. 


2.2.3. Grain-surface network 

For completeness, we also supplement our reaction scheme with 
grain-surface association reactions extracted from the publicly 
available Ohio State University (OSU) network (Garrod et al. 
2008). For those species important in grain-surface chemical re¬ 
action schemes, e.g., the CH 3 0 radical, which are not included 
in Rate 12, we also extract the corresponding gas-phase for¬ 
mation and destruction reactions from the OSU network. The 
grain-surface network has been further updated to include all 
studied routes to water formation under interstellar and circum- 
stellar conditions (jCuppcn et ak] [2010aj Lamberts et al. 2013) . 
The grain-surface reaction rates are calculated assuming the 
Langmuir-Hinshelwood mechanism only, and using the rate- 
equation method as described in Hasegawa et al. ( [1992) . We 
limit the chemically ‘active’ zone to the top two monolayers of 
the ice mantle. We assume the size of the barrier to surface dif¬ 
fusion is 0.3 x the binding energy; in this way, volatile species 
diffuse at a faster rate than strongly bound species. This value 
lies at the optimistic end of the range determined by recent off- 
lattice kinetic Monte Carlo simulations of surface diffusion of 


e.g., Baulch et al. |2005) which are necessary because three- 


CO and C0 2 on crystalline water ice (Karssemeijer & Cuppen 
2014) . This allows the efficient formation of complex organic 
molecules via radical-radical association reactions at > 20 K 


http://www.udfa.net 


http://kinetics.nist.gov/kinetics/index.j sp 
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Table 1 . Star and disk parameters. 


Star 

M* 

(M q ) 

R+ 

(Re) 

T* 

(K) 

M 

(Mo yr 1 ) 

£ioau 
( g cnr 2 ) 

L x 

(erg s^ 1 * ) 

UV Excess? 

M Dwarf 

0.1 

0.7 

3,000 

io- 9 

1.0 

10 3 ° 

Y 

TTauri 

0.5 

2.0 

4,000 

10~ 8 

12 

10 3 ° 

Y 

Herbig Ae 

2.5 

2.0 

10,000 

10~ 8 

56 

10 29 

N 


Notes. The total disk mass will depend on the assumed outer disk radius which in this work, is » 10 AU. Instead, we list the gas mass surface 
density at a radius of 10 AU. 


(see, e.g., fVasy unin & Hcrhst pOl 3[ Walsh et al.|2014 i. For the 


lightest reactants, H and H 2 , we use either the classical diffu¬ 
sion rate or the quantum tunnelling rate depending on which is 


fastest (see, e.g., Tielens & Hagen 1982; 

Hasegawa et al.| 1992). 

For the latter rates, we follow Garrod & 

Yuly (|2011) and adopt 

a rectangular barrier of width 1.5 A. We also include reaction- 
diffusion competition in which the reaction probability is deter¬ 
mined by the relative rates between the barrier-mediated reaction 

and thermal diffusion (see, e.g., Chang et al. 

2007; Garrod & 


Pauly 2011). Although still relatively simplistic, this takes into 
account the increased probability of reaction in the limit where 
the thermal diffusion of the reactants away from a common bind¬ 
ing site is slow compared with the barrier-mediated reaction rate. 


2.2.4. Photodissociation 

One further important process now included is a more robust 
description of the photodissociation rate of N 2 which includes 
the effects of self shielding and mutual shielding by H and H 2 . 
N 2 is important because it is thought to be the main nitrogen¬ 
bearing molecule in interstellar and circumstellar media. Self 
(and mutual) shielding occurs predominantly for those species 
which dissociate via line transitions and occurs when foreground 
material removes photons necessary for dissociation deeper into 
the cloud. In this way, the photodissociation rate of molecules 
which can self shield is reduced relative to the dissociation rates 
for those species which dissociate via the absorption of contin¬ 
uum photons only. H, and CO are famous examples of molecules 


which can self shield (see, e.g., Federman et al.|1979 Glassgold 

et al. 

1985; 

van Dishoeck & Black 

1988 

Lee et al. 1996). H-, also 


et al. (2011): an ‘effective’ shielding column is calculated based 
on the calculated FUV integrated flux relative to that assuming 
no intervening material. The FUV extinction, Tuv. is given by 


t\jv(R,Z) - - In 


Gy\AR, Z) 
GJR,Z) + G C 


( 1 ) 


where G\jy(R,Z) and G*(f?,Z) are the calculated and the geo¬ 
metrically diluted unattenuated stellar FUV integrated fluxes at 
a grid point ( R,Z ), and G ex t is the external unattenuated FUV 
flux. The effective visual extinction is related to the UV extinc¬ 
tion via the empirical relation, A' y (R,Z ) ~ Tyjy(R, Z)/3.02 mag, 
and the ‘effective’ H, column density is calculated using, 
( R,Z ) ~ 0.5 x (1.59 x 10 21 ) X A' W (R,Z) cm~ 2 (Bohlin et al. 
119781 ). For simplicity, the effective shielding column densities 
~f 2 at each point in the disk are estimated by assuming N 2 

io- 5 . 


for N 

has a fixed (rather conservative) fractional abundance, 


of CO and N 2 leading to shielding of CO and N, by foreground 
FU, a process termed ‘mutual’ shielding. 

The inclusion of the self (and mutual) shielding of N 2 is 
now possible due to the work by Li et al. ( 2013a| > in which they 
present parametrised shielding functions calculated using a high- 
resolution model spectrum of N 0 . These shielding functions are 
publicly available for downloacnto use in astrochemical models. 
To use the computed shielding functions which are parametrised 
in temperature and H, H 2 , and N 2 column density, one needs to a 
priori calculate the foreground column densities of each species. 
In protoplanetary disk models this is not trivial because the dis¬ 
sociating radiation can have multiple sources (e.g., stellar pho¬ 
tons and interstellar photons). In addition, the chemistry needs 
to be calculated in series (from the inside outwards and from 
the surface downwards) rather than in parallel which can sig¬ 
nificantly increase computation time. For these reasons, a more 
pragmatic approach is adopted here, similar to that used in Visser 


with respect to H 2 . Finally, the photodissociation rate for N 2 is 
given by, 

k^(R,Z) = *„’(/?, Z)x0 n 2 [ N' n (R,Z),N' n (R,Z),T(R,Z)\ s' 1 , 

( 2 ) 

where k^ 2 (R,Z) = J l Guv(R,Z,A)cr N ^(A)dA is the unshielded 
photodissociation rate for N 2 and f>t\, is the shielding function 
which is a function of H, and N, column density and tempera¬ 
ture, T. 

We adopt a similar method for the computation of the H , and 
CO photodissociation rates using shielding functions calculated 
by|Lee et al.|(|1996|) and|Visser et al.|(|2009 ), respectively. 


2.2.5. X-ray-induced reactions 

We include a set of X-ray-induced reactions which we dupli¬ 
cate from the existing set of cosmic-ray-induced reactions con¬ 
tained in Rate 12 ( McElroy et al.|[2013j ). The reaction rates are 
estimated by scaling the cosmic-ray-induced reaction rate by 
the ratio of the local X-ray and cosmic-ray ionisation rates, i.e., 
k X R ~ U-rX(G<r/4cr)- This is a common assumption in chemical 
models of X-ray irradiated environments (see, e.g., Maloney et 
|al.|1996||Stauber et al.|20 05). The ‘secondary’ X-ray ionisation 
rate, £ X r> is calculated at each grid point in the disk by taking 
into account the local X-ray spectrum and the explicit elemental 
composition of the gas ( |Glassgold et al.|1997[ ). We also include 
a set of ‘primary’ X-ray ionisation reactions (for further details 
see | Walsh et al.|2012| ). 


3 http://faculty.Virginia.edu/ericherb/research.html 
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Fig. 1. FUV radiation fields at a distance of 1 AU from the central star 
(see text for details). 


2.2.6. Initial Abundances 


To generate a set of initial abundances for input into the disk 
model, we run a dark cloud model (7 gas = 7' ( i ust = 10 K, 
n = 10 4 cm~ 3 , and ,4 v = 10 mag). We use the low-metal elemen¬ 
tal abundances from Graedel et al. ( ,1982 1 supplemented with 
updated values for O, C, and N based on diffuse cloud obser¬ 
vations: 3.2 x 10 4 , 1.4 x 10 ~ 4 , and 7.5 x 10 3 relative to total 


hydrogen nuclei density, respectively (Cardelli et al. 1991 1996 


Mey er et al.| [T998 ). In this way, we begin the disk calculations 
with an appreciable ice reservoir on the grain mantle built up 
over the lifetime of the pre-stellar core prior to disk formation. 


In Table[2]we list the abundances (with respect to total H nu¬ 
clei density) of abundant C-, N-, and O-bearing species at times 
of 1.0, 3.2, and 10.0 xlO 5 years. We limit the listed species to 
those which have an abundance >1% that of water ice, the dom¬ 
inant O-bearing species at late times (> 10 5 years). Over the rel¬ 
atively large time steps listed, the trend from atomic to molecular 
gas can be seen, as can the freezeout of volatile species formed 
in the gas phase, such as CO and N-,. Ice species formed in situ 
via hydrogenation of atoms on and within the grain mantle (CH 4 , 
NH 3 , and H.,0) show a general trend of increasing abundance to¬ 
wards late times (~ 10 6 years) as does the abundance of CO-, ice 
which is formed primarily via the reaction between CO and OH. 
The behaviour of H 2 CO and CH 3 OH is more complex: although 
both are formed via the hydrogenation of CO ice, over time, pro¬ 
cessing of the ice by the cosmic-ray-induced UV field causes a 
depletion in CH 3 OH at late times in favour of H 2 CO. 


The median relative ice abundances in dense, quiescent 
clouds measured in absorption against background stars is 


100:31:38:4 for H 2 0:C0:C0 2 :CH 3 0H (Oberg et al. 


20111 . 


These values are also in line with those measured in low-mass 
protostellar envelopes which also include measurements for NH 3 
and CH 4 (H 2 0:NH 3 :CH 4 = 100:5:5, Oberg et al. 2011). We opt 
to use initial abundances at a time of 3.2 xl(r years which cor¬ 
responds to an ice ratio of 100:15:6:1 for H 2 0:C0:C0 2 :CH 3 0H 
and 100:7:22 for H-,0:NH 3 :CH 4 . This is a compromise between 
the set of abundances at early and late times: at the former, the 
CO abundance is low compared with observations, whereas at 
the latter, the CH 3 OH abundance is low compared with obser¬ 
vations. However, for the physical conditions in the inner disk, 
we expect the initial abundances to be important only in the mid- 


Table 2. Abundances with respect to total H nuclei density for the dark 
cloud model and used as initial abundances in the disk model. 


Species 

l.OxlO 5 yrs 

3.2xl0 5 yrs 

l.OxlO 6 yrs 


Gas 

species 


H 

7.7(-05) 

5.2(-05) 

3.61-05) 

h 2 

5.01-01) 

5.0(-01) 

5.01-01) 

He 

9.8(-02) 

9.81-02) 

9.81-02) 

C 

4.4(-05) 

2.2(-06) 

5.11-09) 

N 

5.1(-05) 

2.8(-05) 

5.81-07) 

O 

1.91-04) 

1.0(-04) 

7.51-06) 

ch 4 

1.0(-06) 

8.8(-07) 

6.51-08) 

CO 

5.5(-05) 

7.4(-05) 

1.91-05) 

N 2 

3.2(-06) 

6.0(-06) 

4.2(-06) 

o 2 

8.2(-09) 

1.8(-07) 

5.41-06) 


Ice 

species 


ch 4 

1.7(-05) 

2.51-05) 

2.81-05) 

nh 3 

5.8(-06) 

7.41-06) 

8.41-06) 

h 2 o 

6.1 (-05) 

1.11-04) 

1.81-04) 

CO 

6.91-07) 

1.61-05) 

6.71-05) 

n 2 

2.4(-06) 

1.11-05) 

2.81-05) 

h 2 co 

2.4(-08) 

1.51-07) 

2.31-06) 

ch 3 oh 

1.61-06) 

1.51-06) 

9.31-07) 

co 2 

2.51-06) 

6.31-06) 

1.31-05) 


Notes. a(b) means a x 10* 


plane of the disk where the chemical timescales can be long com¬ 
pared with the disk lifetime. 


3. Results 

3. 1. Disk physical structure 

In Figure [2] we display the physical structure of each disk as 
a function of radius and height (scaled by the radius) for the 
M dwarf disk (left-hand column), T Tauri disk (middle column), 
and Herbig Ae disk (right-hand column). The lower density in 
the atmosphere of the Herbig Ae disk is because the scale height 
of the Herbig disk is smaller than that for the disks around the 
lower-mass stars (77 = c s /Q oc M* 05 , where c s and Q are the 
sound speed and keplerian angular velocity, respectively). The 
surface density of the M dwarf disk is around an order of magni¬ 
tude lower than the other two objects; hence, the lower number 
density of gas throughout. In all three disks, the gas and dust 
temperatures decouple in the disk atmosphere such that the gas 
is significantly hotter than the dust. There are several general 
trends with increasing spectral type: (i) the gas and dust tem¬ 
peratures increase; (ii) the strength of the FUV flux in the disk 
surface increases; (iii) the strength of the X-ray flux in the disk 
surface decreases. For the M dwarf and T Tauri disks, the X-rays 
penetrate deeper into the disk atmosphere than the FUV photons. 
For all three disks, the midplane is effectively shielded from all 
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sources of external radiation, including the central star and the 
interstellar radiation field. The increasing importance of viscous 
heating in the midplane is indicated by a temperature inversion 
below which the temperature begins to increase with depth (see 
the second and third rows of Figure^. 

3.2. C 2 H 2 and HCN 
3.2.1. Chemical structure 

In Figure [3] we display the fractional abundance with respect to 
gas number density of C 2 H 2 and HCN as a function of disk 
radius, R , and disk height divided by radius, Z/R. The dot¬ 
ted and dot-dashed lines represent the t — 1 surface at 3 and 
14 pm, respectively. These are determined using the dust opac¬ 
ity table adopted in the computation of the disk physical struc¬ 
ture. Throughout the remainder of the paper, we adopt the ad 
hoc definition of the disk atmosphere as the material above the 
r( 14 pm) = 1 surface. 

As the stellar effective temperature increases, several trends 
are evident: (i) the disk molecular layer is pushed deeper into 
the disk atmosphere; (ii) the fractional abundances of C 2 Ht and 
HCN decrease in the atmosphere; and (iii) the extent over which 
both species reach a significant abundance increases in the disk 
midplane. C 2 H 2 and HCN are relatively abundant in the molec¬ 
ular layer of the M dwarf disk reaching maximum fractional 
abundances of « 5 x 10 6 and ~ 5 x ICC 5 , respectively. The 
peak C 2 H 2 fractional abundance in the molecular layer of the 
T Tauri disk is ~ 1 x 10 7 , whereas that for the Herbig Ae disk 
is negligible (< 10 11 ). The peak HCN fractional abundance in 
the molecular layer for the T Tauri and Herbig Ae disks is around 
two (~ 10 7 ) and four (~ 10 9 ) orders of magnitude lower than 
that for the M dwarf disk. Thus, the model results suggest that 
the relative molecular complexity in the disk atmosphere de¬ 
creases with increasing stellar effective temperature in line with 
increased photodestruction. 

C 0 H 2 reaches a relatively large fractional abundance (> 10 6 ) 
only in specific regions in the midplane of the T Tauri and Her¬ 
big Ae disks whereas in the M dwarf midplane it is much lower 
(< 10 9 ). In contrast, HCN reaches a relatively high fractional 
abundance (> 10 6 ) over a greater spatial extent when compared 
with C 2 H 2 . An apparent HCN ‘snow line’ moves outwards as the 
disk midplane temperature increases (see Figure [3J. However, 
emission from midplane HCN at near- to mid-IR wavelengths is 
likely obscured by dust in the disk atmosphere as shown by the 
locations of the r = 1 surface at 3 and 14 pm in Figure [3] 


3.2.2. Chemistry of C 2 H 2 and HCN 

A snapshot of the dominant chemical reactions controlling the 
abundance of C-,H 2 and HCN (and related species) in the disk 
atmosphere is given in Figures [4] and [5] respectively. The reac¬ 
tions shown are those which contribute > 10% to the forma¬ 
tion and destruction rates at the position of peak fractional abun¬ 
dance in the atmosphere at a radius of 1 AU. The networks are 
similar to those presented in Agundez et al. ( 2008[ > and Bast et 
[aT7| ( |2013| ) except that we also include the dominant destruction 
mechanisms. 

Free carbon and nitrogen (necessary for incorporation into 
molecules and radicals such as CH, NH, CN, and C 2 ) are re¬ 
leased from the main gas-phase reservoirs (CO and N,) via pho¬ 
todissociation for the Herbig Ae disk. For the two cooler disks, 
both X-ray-induced dissociation and reactions with He + (which 
itself is produced by direct X-ray ionisation) are dominant. Many 
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reactions involving He have measured rate coefficients (Adams 
|& Smith| 1979j |Anicich et al.|1977) . 

Whether or not a molecule survives in the disk atmosphere 
with an appreciable abundance requires a delicate balance be¬ 
tween formation (via neutral-neutral chemistry or ion-molecule 
chemistry) and destruction (via photodissociation or X-ray- 
induced dissociation). The abundance and distribution of HCN 
in the atmosphere of all three disks is primarily controlled by 
CN (|Baulch et al.119941, a 
al. 20061). 


formation via H , 


photodissociation (|van Dishoeck et 


and destruction via 
2006). In the Herbig Ae 


disk, destruction via reaction with atomic hydrogen also plays a 
minor role. This reaction has a large reaction barrier (12,500 K, 
Tsang & Herron] [T991) and is only significant in very hot gas 
(> 1000 K). CN has numerous formation routes via neutral neu¬ 
tral reactions: N + C 2 ((Smith et al.|2004), NH + C (jBrownsword 
|et al.|[T996| ), and NO + C ( |Chastamg et al.||2000[) wh ere NO is 
formed via the reaction, N + OH (Wakelam et al. 2012). Only the 
formation of NH (via N + H 2 ) possesses a substantial reaction 
barrier (18,095 K, |McElroy et al.|2013) . 

For the M dwarf and T Tauri disks, photodissociation at 
Lyman-a is more significant than that by the FUV continuum 
background because around 70 - 80% of the FUV flux is con¬ 
tained in the Lyman-a- line (see also Fogel et al.|20TT I. This per¬ 
centage is in line with that determined towards observations of 
M dwarf and classical T Tauri stars (see, e.g., France et al.|2013| 
2014|l. HCN is treated as though it photodissociates via line tran¬ 


sitions (|Lee|1994] see, e.g.,) and it also has a non-negligible pho 


van 


todissociation cross section at 1216 A (<x =3 x 10 ' cm 2 
[Dishoeck et al.|2006| ). 

HCN is relatively more abundant in the M Dwarf disk due to 
the weaker FUV flux leading to decreased destruction via pho¬ 
todissociation. Ion-molecule reactions also contribute to HCN 
formation in the M dwarf disk (via HCNH + + e ) and destruc¬ 
tion (via HCN + X + , where X + is H + , C + , H}, H 3 0 + , and He + ) as 
shown in Figure 5] These reactions all have measured rate coeffi¬ 


cients (|Huntress 1977| Clary et al.| 1985 1 Anicichetal.|1993 Se- 
maniak et al.|2001) . X-ray-induced photodissociation also con¬ 

tributes to the destruction of HCN at the level of 15- 20%. The 
ion-molecule formation route is triggered by the formation of 
CN + and HCN + via the reactions, N + CH + , NH + C + , and 
N + CH} (|Viggiano et al.|1980||Prasad & Huntress|1980|). HCN + 

+ H , and 


HCN 


H, 


and HCNff are then formed via CN 4 
( |Raksit et al.|1984[|Huntress|1977| ). 

In all cases, midplane gas-phase HCN is synthesised via 
neutral-neutral chemistry: it is not related to the desorption of 
HCN ice as is the case for traditional snow lines. The abun¬ 
dance is mediated by format ion via H 2 and CN and destruction 
via collisional dissociation (Baulch et al.|1994||Tsang & Herron 
|1991j >. This formation route has a reaction barrier of 820 K and 
thus requires warm temperatures for activation (> 200 K). These 
temperatures are surpassed in the midplane of each disk due to 
the inclusion of heating via viscous dissipation (for details see 
|Nomura & Millar||20051 . and references therein). Viscous heat¬ 
ing dominates over stellar heating in the midplane within radii 
of ~ 0.40, 2.5, and 1.1 AU for the M dwarf, T Tauri, and Her¬ 
big Ae disks, respectively. The Herbig Ae transition radius is less 
than that for the T Tauri disk because the former has significantly 
stronger stellar heating. 

In Figured] a snapshot of the chemistry of C 2 H 0 is presented. 
In all three disks, the formation and destruction of C-,H 2 in the at¬ 
mosphere is dominated by the neutral-neutral reaction H, + C 2 H 
(jLaufer & Fahr|2004 1 and photodissociation (van Dishoeck et al. 
2006|), respectively. C 2 H 2 is also preferentially photodissociated 
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Fig. 2. Number density (top row), gas temperature (second row), dust temperature (middle row), far-UV flux (forth row) and X-ray flux (bottom 
row) as a function of disk radius and disk height (scaled by the radius) for the M dwarf disk (left column), the T Tauri disk (middle column) and 
the Herbig Ae disk (right column). 


van 


at Lyman-o' wavelengths in both cases (<x > 4 x 10 17 cm 2 , 
Dishoeck et al.|[2006|>, for similar reasons as discussed above 


for HCN. In the M dwarf disk, C,H, 
actions with C, H{ 


also destroyed via re 


and He (Kim & Huntress 1975 Chastaing 


et al.|199~9j|Laufer & Fahr|2004| i, as well as via X-ray-induced 

photodissociation (at the level of » 15 - 20%). There are also 
formation routes via ion-molecule chemistry (C 2 H J + e ) which 
contribute at the level of a few percent and are barrierless. 


C 2 H 2 is significantly less abundant in the atmosphere of the 
Herbig Ae disk than in the other two disks. In the absence of ef¬ 
ficient ion-molecule pathways (triggered by the formation of C + 
via CO + He + ), an important first step in the formation of carbon- 
chain molecules is the formation of CH via C + H 2 which has a 
large reaction barrier (11,700 K). CH can then react barrierlessly 
with atomic C to give C 2 . In the atmospheres of protoplanetary 
disks, the gas temperature is controlled by the strength of the UV 


Article number, page 7 of|28| 













































































































A&A proofs: manuscript no. inner_disk_chemistry_accepted 



0.1 


1 

R (All) 


10 0.1 


1 

R (AU) 


10 0.1 


1 

R (AU) 


10 


Fig. 3. Fractional abundance relative to total gas number density of C,H, (top row) and HCN (bottom row) for the M dwarf disk (left-hand 
column), T Tauri disk (middle column), and Herbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density 
(integrated from the surface downwards) at which r « 1 at 3 pm and 14 pm, respectively. 


field. Carbon-chain growth is impeded in the atmosphere of the 
Herbig Ae disk because the increased photodissociation counter¬ 
acts the temperature-activated gas-phase chemistry. In Figure [4} 
we split the carbon chemistry into that dominated by UV radia¬ 
tion (and a higher gas temperature) highlighted in red on the left- 
hand side, and that dominated by X-ray radiation (and a lower 
gas temperature) highlighted in blue on the right-hand side. The 
left-hand side represents the chemistry more dominant in Her¬ 
big Ae disks and the right-hand side represents the chemistry 
more important in M dwarf disks. The chemistry tends towards 
the middle of this reaction scheme in environments where the 
UV radiation is too strong for the survival of species other than 
CO, C, and C + . 

In summary, the gas-phase chemistry depends, not only on 
the strength of the FUV radiation (which controls the gas tem¬ 
perature), but also on the adopted ionisation sources and spectra 
and corresponding rates propagated throughout the disk. This 
is especially true for the M dwarf disk in which ion-molecule 
chemistry also plays a role in the formation and destruction of 
HCN and C 2 H 0 as shown in Figures |4] and [5] The liberation of 
free carbon and nitrogen is the principle underlying reason for 
the importance of X-ray chemistry. 

3.2.3. Column densities 

A higher C 2 H 2 and HCN fractional abundance is seen in the 
M dwarf disk atmosphere compared with the others; however, 
this does not necessarily translate into an observable column 
density, especially given the more tenuous nature of the M dwarf 
disk and the fact that the molecules peak in fractional abundance 
higher in the disk atmosphere (where the density is also lower). 
Figure [6] displays the vertically integrated column densities of 
both species as a function of disk radius over the entire ver¬ 
tical extent (left-hand column) and down to the r = 1 surface 
at 14 pm (right-hand column). There is a general trend that the 
column density peaks at inner radii, then falls off sharply at a 
particular radius which moves outwards with increasing stellar 



Fig. 4. Dominant reactions controlling the abundance of CkH 2 and re¬ 
lated species at the position of peak fractional abundance in the disk 
atmospheres at a radius of 1 AU. The solid lines indicate gas-phase 
two-body reactions and the dotted lines indicate UV or X-ray-induced 
reactions. The width of the arrows for the gas-phase reactions shows 
the magnitude of the activation barrier. The left-hand side, highlighted 
in red, is the chemistry most dominant in UV-dominated gas with a 
higher gas temperature (e.g., in Herbig Ae disks), whereas that high¬ 
lighted in blue on the right-hand side is that which dominates in X-ray 
irradiated gas with a lower gas temperature (e.g., in X-ray irradiated 
M dwarf disks). 
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Fig. 5. Dominant reactions controlling the abundance of HCN and re¬ 
lated species at the position of peak fractional abundance in the disk 
atmospheres at a radius of 1 AU. The solid lines indicate gas-phase 
two-body reactions and the dotted lines indicate UV or X-ray-induced 
reactions. The width of the arrows for the gas-phase reactions shows the 
magnitude of the activation barrier. 


spectral type. This fall off occurs at a smaller radius for C 0 H 2 
than for HCN. This behaviour is also seen in the T Tauri model 
presented in |Agundez et al.||2008) l. 

For C-,Ht. the T Tauri and Herbig Ae disks achieve a similar 
peak column density ~ 10 20 cm~ 2 at radii of ~ 0.2 and 0.3 AU 
respectively. The column density then falls with increasing ra¬ 
dius to < 10 17 cm~ 2 beyond 0.5 AU for the T Tauri star and 
beyond 1 AU for the Herbig Ae star. The total column density 
for the M dwarf disk lies orders of magnitude lower, reaching a 
peak value of ~ 10 17 cm -2 at 0.1 AU and falling to < 10 15 cm -2 
beyond 1 AU. Comparing with the column densities calculated 
down to the r(14 fin l) = 1 surface, for the two warmer disks, 
the values are lower by between three and four orders of magni¬ 
tude because the total column density is dominated by midplane 
C 2 H 2 . In contrast, the values for the M dwarf disk remain com¬ 
parable beyond ~ 1 AU because in this case, the total column is 
dominated by atmospheric C 2 H 0 . The M dwarf disk has a larger 
column density in the atmosphere than the T Tauri disk beyond 
0.2 AU. Although C 2 H, does not appear to be abundant in the 
atmosphere of the Herbig Ae disk (according to Figure [3j, there 
is still a significant column density: this is because there is a thin 
layer (only one to two grid cells wide) of relatively abundant 
C 0 H 2 which overlaps with the r = 1 surfaces. This is likely a 
feature of the grid resolution of our model and thus is a numeri¬ 
cal artefact. 

For HCN, the T Tauri and Herbig Ae disks reach a peak col¬ 
umn density of ~ 10 21 cm 2 at similar radii to those for C 2 H-,. 
The radial behaviour of the column density then follows the spa¬ 
tial extent of HCN in the disk midplane (see Figure[3]). Similar to 
C 0 H 2 , the peak value for the M dwarf disk is lower (~ 10 20 cm~ 2 
at 0.1 AU), and the column density remains constant beyond 
~ 0.8 AU at a few x 10 15 cm -2 . For the corresponding values 
down to the t(14 fim) = 1 surface, the M dwarf disk has a higher 
column density than the T Tauri and Herbig Ae disks beyond 
0.6 AU and 4 AU, respectively. Although the fractional abun¬ 
dance of HCN is lower in the Herbig Ae disk, the molecular 
layer is located deeper in the disk atmosphere where the den¬ 
sity is higher leading to the apparently large calculated column 
density of HCN (for "R <4 AU). 


3.2.4. Comparison with observed trends 

The preferred method for comparing model results with obser¬ 
vations is the simulation of the line emission; however, this is a 
non-trivial matter involving careful consideration of collisional 
and radiative excitation in the IR (see, e.g., [Pontoppidan et al. 
2009| ~ ’ " 

2015 


Meijerin k et al.|[2009l |Thi et alT] |20 13 [ jBruderer et al. 
Moreover, the dust properties and size distribution (in¬ 
cluding, e.g., grain growth) become crucial for the calculation of 
the emitted spectrum (see, e.g., |Meijerink et ak|2009| . This is 
beyond the scope of the work presented here which is focussed 
on the chemistry; however, this is planned future work. Here, we 
compare the calculated column densities (and ratios) with those 
derived from observations to investigate if the chemical models 
are at least able to reproduce the observed abundances and re¬ 
lated trends. 

The C 2 H 2 /HCN column density ratio derived for disks 


around cool stars ranges from 0.87 to 4.3 (see Table 6 in Pas 


cucci et al. 20131. In contrast, the C 2 H,/HCN column density 
ratios derived for T Tauri disks lies between 0.006 and 0.43 (see 


Table 4 in Carr & Najita 2011 1 . Salyk et al. (2011 1 derive a range 
from 0.13 to 20 for their sample of T Tauri stars. In Figure [7] 
the model ratios in the disk atmosphere are plotted as a function 
of radius for the T Tauri disk (red dashed lines) and M dwarf 
disk (gold dotted lines), overlaid with the observed range. The 
fine lines are the equivalent ratios for the column densities inte¬ 
grated down to the r = 1 surface at 3 fim. The ratio is flat for the 
M dwarf disk (~ 0.1) and lies roughly one order of magnitude 
lower than the observed range (gold hatched zone). On the other 
hand, the ratio for the T Tauri disk increases with radius reaching 
a peak of ~ 4. Also, the T Tauri ratios beyond 0.7 AU lie well 
within the wide range of observed values (red shaded zone). 

For HCN, Carr & Najita ( |2011| derive best-fit column 
densities for their sample of T Tauri disks between 1.8 and 


6.5 x 10 ’ cm , however, assuming extremes in the optical 
depth of the HCN emission expands this to between 0.2 and 
31 x 10 16 cm~ 2 . F or C 2 H 2 , t he derived range is 0.02 to 1.6 x 
10 16 cm -2 . Salyk et al. (201 lj) derive column densities between 

\ Q V 1 H r , rr - t ~~2 nrwj V\otu;oan f) onrl 1 A V 1 A^ 


2 and between 0.05 and 1.0 x 10 16 cnD 2 
for HCN and C 2 H 2 respectively. The observed column density 
ranges are indicated by the solid colour regions in Figure [6] The 
T Tauri model results give good agreement with both the abso¬ 
lute column densities and the ratio of C 2 fF/HCN in the vicinity 
of the radius within which the observed emission originates. 

In the inner region of the disk (< 2 AU), the models qualita¬ 
tively reproduce the trend that the C 2 H 2 /HCN ratio is higher in 
M dwarf disks than in T Tauri disks; however, the ratio predicted 
in the model lies lower than that observed. Hence, the M dwarf 
model is either overpredicting the HCN column density or un¬ 
derpredicting the C 2 Ht column density. Pascucci et al. (2013) 


derive column densities ranging from 4.5 to 50.1 x 10 16 cm 2 

2 for HCN. The range 


for CtH 2 and from 1.1 to 18.1 x 10 16 cm 


of observed column densities are indicted by the gold hatched 
regions in Figure |6j The model values for C 0 H 2 are more than 
one order of magnitude lower whereas the values for HCN lie 
within a factor of a few of the observed column densities. How¬ 
ever, the absolute values will depend to a degree on the adopted 
model parameters, such as disk surface density. It is also pos¬ 
sible that the excitation mechanisms for C 2 H 2 and HCN differ 
in cool stars relative to T Tauri stars which means that the rel¬ 
ative line emission no longer traces the relative abundances nor 
column densities. The excitation of HCN in protoplanetary disks 
was investigated in detail by Bruderer et al. ( [2015) 1 who con¬ 
clude that HCN abundances derived assuming LTE should differ 
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Fig. 6. Column density of CiHo (top row) and HCN (bottom row) as a function of radius for each disk model for the entire vertical extent of the 
disk (left-hand column) and down to the r = 1 surface at 14 pm (right-hand column). The solid coloured and hatched regions indicate the range of 
observed values for the T Tauri (Carr & Najita 2011, Salyk et al. 2011) and M dwarf (Pascucci et al.|20f3| disks, respectively. 


by no more than a factor of three from those derived assuming 
non-LTE. They also conclude that C 2 H n will behave similarly to 
HCN because of the presence of similar infrared bands through 
which the excitation can be pumped. We plan to test this via sim¬ 
ulations of the molecular emission in future work. 


Pascucci et al. (12013}) speculate that the self-shielding of N 2 
may play a role in determining the C 2 H 2 /HCN ratio: this can 
potentially lock up a greater fraction of atomic nitrogen thereby 
impeding the production of other nitrogen-bearing species, such 
as HCN. Here, we have also shown that X-ray-induced chemistry 
may play an important role in both the production and destruc¬ 
tion of molecules in the atmosphere of disks around cool stars. 
We further investigate the role of N 2 se lf-shielding and X-ray- 
induced chemistry in Sects. |4~T| and |4.2| respectively. 


3.3. OH and H 2 0 

3.3.1. Chemical structure 


Another interesting trend seen in the IR data is the lack of water 
detections in disks around hotter stars. In Figure[8]we display the 
fractional abundance of OH (top row) and H 2 0 (bottom row) as 
a function of radius ( R ) and height divided by the radius (Z/R). 
OH is more extended and resides in a layer slightly higher in 
the disk than HCN, C 0 H 2 , and H 2 0, in line with the hypothe¬ 
sis that the OH chemistry is driven by photodissociation. The 
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peak fractional abundance of OH is 8, 5, and 4 x 10 for the 
M dwarf, T Tauri, and Herbig Ae disks, respectively. The corre¬ 
sponding values for H.,0 are 5, 5, and 4 x 10 4 , again showing 
a general (albeit very shallow) decline in molecular complexity 
with increasing spectral type. The fractional abundance of OH is 
negligible (< 10 11 ) in the midplane of all three disks. The ex¬ 
tent over which H^O is abundant in the disk midplane is related 
to the thermal desorption of H 2 0 ice (we assume a binding en¬ 
ergy of 5570 K, Fraser et al.|2001|). The H 2 0 ‘snow line’ shifts 
outwards in radius with increasing stellar spectral type with po¬ 
sitions at 0.35, 1.5, and 6.1 AU for the M dwarf, T Tauri, and 
Herbig Ae disk, respectively. Beyond these radii, H.,0 is frozen 
out on dust grains in the midplane. 


3.3.2. Chemistry of OH and H 2 0 

Figure [9] shows a snapshot of the dominant reactions contribut¬ 
ing to the formation and destruction of OH and H 2 0 (and related 
species) at the position of peak abundance in the disk atmosphere 
at a radius of 1 AU. Gas-phase H 2 0 is predominantly produced 
in the disk atmosphere via the neutral-neutral reaction, H 2 + OH 


(using the rate coefficient from Oldenborg et al. 1992) . In the 
cooler M dwarf disk, thermal desorption of water ice dominates 
at small radii (~ 1 AU). Destruction is typically via photodissoci¬ 
ation: for the M dwarf and T Tauri disks, this is again dominated 
by Eyman-cr photons (van Dishoeck et al. 2006[l. In the M dwarf 
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Fig. 8. Fractional abundance relative to total gas number density of OH (top row) and H,0 (bottom row) for the M dwarf disk (left-hand col¬ 
umn), T Tauri disk (middle column), and Herbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density 
(integrated from the surface downwards) at which r » 1 at 3 /im and 14 pm, respectively. 



Fig. 7. Ratio of C-H, and HCN column densities in the disk atmosphere 
as a function of radius for the M dwarf disk (gold dotted lines) and 
T Tauri disk (red dashed lines). The solid coloured and hatched regions 
correspond to the range of observed values for T Tauri (Carr & Najita 
|201 1| |Salyk et al.|2011| ) and M dwarf (Pascucci et al.|2013|) disks, re¬ 
spectively. The thick and fine lines represent the ratios down to the r = 1 
surface at 14 and 3 pm. 


disk, there are additional destruction routes via ion-molecul e re¬ 
actions with H + and H 2 (Kim et al. 1974 Smith et al. 1992 1 and 
X-ray-induced photodissociation contributes at the level of 12 - 
15% (similar to that found for C n H 2 and HCN). Similarly, the 
gas-phase abundance of OH is mediated by production via the 
photodissociation of H z O and the reaction between H 2 and O 


Baul ch et al.| 1 992), with destruction via photodissoc iation (|van 
Dishoeck et al. 2006) and reactions with H 2 and H (Tsang et al. 


986 |Qldenborg et al. I992jl. As in the case for H 2 0, dissoci- 


ation by Lyman-a dominates over that by the background FUV 
continuum. This can be summarised in the following, rather suc¬ 


cinct, reaction scheme. 


O 


H./iv 


R 2 freezeout 

OH H 9 0 


H, hv 


desorption 


^2®ice 


(3) 


(see also Bethell & Bergin||2009 Woitke et al.] |2009| and Fig¬ 
ure^ i. Thus, the water chemistry in the disk atmosphere is rather 
simple with the relative abundances of O, OH, and H 2 0 con¬ 
trolled primarily by the relative abundances of H and H 2 , the 
gas temperature (which needs to be sufficiently high to activate 
the neutral-neutral chemistry) and the strength of the FUV ra¬ 
diation field (necessary for photodissociation). We find that re¬ 
actions with vibrationally excited or ‘hot’ H, (see Sect. 2.2.1 1 


increase in importance towards the upper atmosphere. However, 
these reactions do not contribute significantly at the position of 
peak abundance and deeper and hence, do not influence the total 
column density of n either OH nor H 2 0. For a detailed review of 
water chemistry see van Dishoeck et al. (2013)>. 


3.3.3. Column densities 

Figure[T0]shows the vertically integrated column densities of OH 
(top row) and H,0 (bottom row) as a function of radius over the 
entire vertical extent of the disk (left-hand columns) and down 
to the r = 1 surface at 14 pm (right-hand column). The column 
densities show less structure compared with those for C-,H 2 and 
HCN. Because the majority of OH is in the disk atmosphere, 
the column densities for both cases are similar and range be¬ 
tween a few times 10 15 to 10 17 cm~ 2 . The OH column density 
also does not vary greatly with radius. The Herbig Ae disk gen¬ 
erally has the largest column density and the M dwarf has the 
lowest, although the values for the T Tauri and M dwarf disks 
are similar beyond as 1.6 AU. All three disks show very high 
total column densities of H 0 O (> 10 21 cm 2 ) within each respec¬ 
tive snow line. The H 2 0 column densities also show a similar 
trend to that for C 7 H 2 and HCN. The column densities start high 
in the inner disk then fall off sharply at a distinct radius for each 
disk (which moves outwards with increasing spectral type). This 
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(15 - 20 AU) than that expected at shorter wavelengths (see also 
the analysis in |Fedele et al.|[2013| >. 

The Herbig Ae model appears to overproduce gas-phase wa¬ 
ter in the inner disk atmosphere relative to observations; how¬ 
ever, the derived upper limits are dependent on the chosen emit¬ 
ting radius and gas temperature. Figure m shows the ratio of 
OH to HtO column density as a function of disk radius for the 
T Tauri model (blue dashed lines) and the Herbig Ae model 
(green dot-dashed lines). The fine lines represent the equivalent 
ratios for the columns integrated down to the r(3 pm) = 1 sur¬ 
face. The model T Tauri ratios lie within the observed range for 
radii less than « 1 AU. The ratio for the Herbig Ae model has a 
value ~ 10 5 - 10 2 within ~ 8 AU, beyond which it increases 
and quickly tends to 1 at larger radii. The analysis of the non¬ 
detections of OH and H,0 in Salyk et al. (2011[> suggest a lower 


Fig. 9. Dominant reactions controlling the abundance of OH and H,0 
and related species at the position of peak fractional abundance in the 
disk atmospheres at a radius of 1 AU. The solid lines indicate gas-phase 
two-body reactions and the dotted lines indicate UV or X-ray-induced 
reactions. The width of the arrows for the gas-phase reactions shows the 
magnitude of the activation barrier. 


limit to the 0H/H ,0 ratio of a few xlO 2 assuming an emitting 


radius which ranges from a few AU to several 10’s of AU. Fedele 
|etaL| ( |20n) suggest a ratio of > 1-25 within an emitting radius 
as far out as 30 AU based on near-IR data. The model values are 
on the cusp of the Salyk et al. ( 201 1|> lower li mit in the very in¬ 
ner region and tend towards the Fedele et ah] (|2011|l lower limit 
at 10 AU. 


radius generally lies beyond that for C 0 H 2 and HCN. Similar to 
the case for HCN, gas-phase HiO in the disk midplane is likely 
obscured by dust at near- to mid-IR wavelengths. The bottom 
right-hand panel shows that the column density of ‘visible’ H 2 0 
is on the order of a few times 10 19 cm -2 within the snow line 
for the T Tauri and Herbig Ae disks. The column density drops 
to ~ 10 16 cm 2 beyond 1.5 AU for the T Tauri disk. The value 
for the M dwarf model remains constant over most of the radial 
extent of the disk at ~ 10 17 cm -2 . 


18 


cm 

17 


The lack of hot water emission from the innermost regions 
of disks around Herbig Ae stars remains a puzzle, especially 
given the hypothesis that water can shield itself from pho¬ 
todissociation by the stellar radiation field at column densities, 
A(H 2 0) > 2 x 10 17 cm -2 (assuming a photodissociation cross 
section, cth,o ~ 5 x 10 18 cm 2 ; Bethell & Bergin 2009 Du 


3.3.4. Comparison with observed trends 

Here, we compare the column densities calculated by the mod¬ 
els with those derived from the observations to test whether the 
chemical calculations are able to reproduce the observed abun¬ 
dances and trends. For their sample of classical T Tauri stars 
in which both OH and H 2 0 were detected, Salyk et al. ( 2011[ > 
derive OH column densities from 0.04 to 6 X 10 1 (1 cm 2 . The 
H^O column densities range from 4 x 10 17 to 1.8 x 10 


to 


Carr & Najita (2011) derive a wider range of values (4 x 10 
7.9 x lO 20 cm -2 ). The observed ranges are indicated by the solid 
colour shaded regions in Figure [TO] Despite the model T Tauri 
disk not being representative of any particular source, there is 
significant overlap between the calculated column densities in 
the disk atmosphere and observed values for both species. 

For the Herbig Ae disks, only OH has been robustly detected 
at IR wavelengths ( [Fedele et al.||2012] with one exception dis¬ 
cussed below). The range of observed column densities (derived 
from VLT/CRIRES data) lie between 1.3 and 20 x 10 16 cm 2 and 
the model values lie within a factor of a few of this range (indi¬ 
cated by the striped region in Figure [T0| [Fedele et al.|2011) >. The 
maximum model H 2 0 column density in the disk atmosphere 
is more than an order of magnitude larger than the maximum 
upper limits derived by both Saly k et ak] ( |2011[ > and Fedele et 
al. (2011 1 . OH and H 2 0 have been detected in the disk around 
HD 163296 in HerschelfPA.CS data at far-IR wavelengths with a 
relative column density of OH/H 2 0 as 1 and an excitation tem¬ 


perature ranging from 200 - 500 K (Fedele et ak |2012] >. How¬ 
ever, this emission originates further out in the disk atmosphere 


& Bergin|2014 Adamkovics et al.|2014 i. The column densities 

predicted here for the Herbig Ae disk are in line with those by 
other work (see, e.g., Woitke et ak 2009) >. One possible expla¬ 
nation is that turbulent mixing within the planet-forming region 
can help sequester water in the midplane where, if the temper¬ 
ature is sufficiently low, it can become trapped as ice on dust 
grains. This would require that as the dust grain grow, either via 
ice mantle growth or via coagulation, they become decoupled 
from the gas and remain in the shielded (and cold) midplane (see, 
e.g., [Stevenson & Lunine [T988] >. There are two caveats to this 
theory: the first is that the temperature of the midplane would 
need to be below « 150 K within ~ 10 AU. The midplane tem¬ 
perature in our Herbig Ae model is generally too high for water 
to reside as ice on dust grains except beyond ~ 6 AU (see Fig¬ 
ure]^; however, the exact temperature profile of the disk is sen¬ 
sitive to numerous factors including, for example, the adopted 
disk surface density, dust-grain size distribution, and degree of 
flaring (see later). It also remains to be demonstrated whether 
large-scale radial mixing between the warm inner midplane and 
cold outer midplane is a viable mechanism in disks around hotter 
stars. The second caveat is that one would expect the abundance 
of OH to also be affected by vertical mixing in the atmosphere 
since OH and HiO are chemically coupled. For a given radiation 
field, the rate coefficients for the primary routes to formation 
and destruction are known; hence, there may be additional de¬ 
struction routes for gas-phase water, not yet included in the net¬ 
works. One potential explanation which remains to be explored 
in protoplanetary disks is the effect of rotationally excited OH 
which can be produced via photodissociation of H.,0 at Lyman- 
a wavelengths (see, e.g., Fillion et a h[2001 1 . The rate coefficient 
for the OH(v, j) + H reaction can be enhanced by several or¬ 
ders of magnitude relative to that for ground state OH (Li et ak 
2013b[ >. These reactions may be important for shifting the ratio 
of O/OH/H-,0 in the disk atmosphere (see Equation]^. 
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Fig. 10. Column density of OH (top row) and H 2 0 (bottom row) as a function of radius for each disk model for the entire vertical extent of the 
disk (left-hand column) and down to the r = 1 surface at 14 pm (right-hand column). The solid coloured and striped regions indicate the range of 
observed values for the T Tauri (|Carr & Najita|201 f||Salyk et al.|201 1|) and Herbig Ae jFedele et al.|201 l[|Salyk et al.|201 1|> disks, respectively. 


|Pontoppidan et al. (2010 1 also extensively discuss several hy¬ 
potheses for the lack of hot H 2 0 in Herbig Ae disks, including (i) 
an intrinsically lower abundance by an, as yet, unknown physi¬ 
cal or chemical mechanism; (ii) veiling of the molecular features 
by the strong mid-IR background (if the water line luminos¬ 
ity is a weaker function of stellar spectral type than the contin¬ 
uum); and (iii) a well-mixed disk atmosphere with the canonical 
gas-to-dust mass ratio of 100. The results presented here show 
that for a well-mixed atmosphere, the observable column den¬ 
sity of water vapour remains high. However, a second outcome 
of grain growth in protoplanetary disks (in addition to that dis¬ 
cussed above) is the greater penetration of FUV radiation which 
can push the molecular layer deeper into the disk atmosphere 
(see, e.g., [Aikawa & Nomura [2006| l; hence, a higher fraction of 
the gas-phase water may be ‘hidden’ from view. This remains 
to be confirmed specifically for the inner regions of disks. An¬ 
other factor to consider, and mentioned previously, is the disk 
gas and dust structure. Disks around Herbig Ae/Be stars have 
been classified into groups based on the shapes of their SED at 
mid-IR wavelengths (Meeus et al. |2001| >. Group I disks are pos¬ 
tulated to have a flared structure which allows the disk to capture 
more FUV photons which increases the gas and dust temper¬ 
ature throughout the disk, whereas group II disks are ‘flatter’, 
capture less FUV and are thus much colder. Dust grain growth 
and settling has been postulated as the reason behind the appar¬ 
ent dichotomy of Herbig Ae/Be disks (see, e.g., |Dullemond &| 


Dominik|2004j > and it has also been suggested that group I disks 


may be transitional in nature, i.e., they have evidence (usually 
confirmed by spatially resolved imaging) of a significant gap in 
the inner disk (see, e.g., Grady et al.|2015) . The disk model we 
have used assumes a flared disk in hydrostatic equilibrium with¬ 
out a gap: whether significant grain growth and the presence of 
an inner gap affects the abundance and distribution of gas-phase 
water remains to be explored. 


4. Discussion 

4.1. On the importance of N z seif shielding 

To quantify the importance of N 2 shielding in the disk atmo¬ 
sphere, in Figure [12] we show the fractional abundance of N, 
N 2 , and HCN as a function of Z/R at R = 1 AU. An equivalent 
figure for R — 10 AU is shown in the Appendix (Figure [AT) . 
Results from three different models are presented: (i) the fidu¬ 
cial model (black dashed lines, with N 2 shielding and X-rays 
included); (ii) a model without N, shielding (purple lines); and 
(iii) a model without X-ray-induced chemistry (green lines). The 
inclusion of shielding has an affect on the relative abundances of 
N and N 2 in a narrow region of the disk only at both 1 and 10 AU 
for all disks. The ratio of N/N 2 is generally more affected as the 
central stellar effective temperature increases. Although the frac¬ 
tional abundance of atomic nitrogen can vary by more than one 
order of magnitude, this translates to a difference on the order 
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Fig. 11. Ratio of OH and H.,0 column densities in the disk atmosphere 
as a function of radius for the T Tauri disk (blue dashed lines) and Her- 
big Ae disk (green dot-dashed lines). The solid coloured region corre¬ 
sponds to the range of observed values for T Tauri disks ( |Carr & Najita 
1 201 I [|Salyk et al.|20 1 1 J and the green diamond and arrow indicates the 
ratio range observed for HD 163296 ( |Fedele et al.|2012) . The thick and 
fine lines represent the ratios down to the r = 1 surface at 14 and 3 pm. 


of a factor of a few only for the HCN abundance. The results 
demonstrate that N 0 shielding alone is not able to account for 
the change in C 2 H 2 /HCN column densities and line flux ratios 
seen from M dwarf to T Tauri stars. 

The above conclusion holds for models in which the dust and 
gas are assumed to be well mixed. If a significant fraction of the 
dust has grown and settled to the midplane (see, e.g., |Dominik etj 
al.|2005 1, this can lead to a relatively dust-poor disk atmosphere 
and allow greater penetration of FUV radiation. In that case, the 
importance of molecular (or self) shielding increases relative to 
dust shielding ( jVisser et al.|2009 |Li et al.|2013a i. Whether the 
same conclusion holds for M dwarf disks with advanced grain 
growth and settling remains to be confirmed. 


4.2. On the importance of X-ray-induced chemistry 

In Figure [l2| (and Figure |AJ~| i models with and without X-ray- 
induced chemistry are also shown (black dashed lines versus 
green solid lines). X-ray-induced chemistry is significantly more 
important for the M dwarf and T Tauri disks than for the Her- 
big Ae disk. In both cases and at both radii, the abundance of 
N 2 is significantly increased in the disk atmosphere when X- 
ray chemistry is neglected which leads to a decrease in the ratio 
of N/N 2 . The abundance of HCN is also significantly perturbed 
by the exclusion of X-ray chemistry. The abundance of HCN is 
increased in the disk atmosphere despite the reduction in avail¬ 
able free nitrogen indicating that X-rays are important for HCN 
destruction even at heights where the FUV field is strong. The 
higher penetration depth of X-rays versus FUV photons allows 
an increase in N/N 2 deeper into each disk and hence leads to 
an increase in HCN in the molecular layer (see also Aikawa & 
Herh stp999] ). For the T Tauri disk at 1 AU, X-rays also help to 
destroy HCN deeper down towards the disk midplane. 

The results show that the inclusion or exclusion of X-ray- 
induced chemistry can have a profound affect on the position 
and the value of peak fractional abundance of molecules such 
as HCN. To investigate whether carbon-bearing species are sim¬ 


ilarly affected, in Figure [13) the fractional abundance of G 2 H 2 
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is shown as a function of Z/R at a radius of 1 AU for both 
the M dwarf and T Tauri disks with and without X-ray-induced 
chemistry (black dashed lines and orange solid lines, respec¬ 
tively). The results show that X-rays are important in the disk 
molecular layer for releasing free carbon into the gas phase for 
incorporation into species such as CtH 2 . The inclusion of X- 
rays in the M dwarf disk increases the peak abundance of C n H 2 
by more than two orders of magnitude and increases the extent 
over which C 2 H, is relatively abundant > 10 9 with respect to 
gas number density. The results for the T Tauri disk are even 
more extreme with a three orders of magnitude increase in the 
peak fractional abundance from « 10 11 to > 10 8 when X-ray- 
induced chemistry is included. 

As discussed in Sect. |3.2| the efficacy of the X-ray-induced 
chemistry is because of the generation of He + which in turn 
reacts with those molecules robust to photodissociation, CO 
and N 2 . This creates a steady supply of free and reactive 
atomic atoms and ions for incorporation into other molecules 
via more traditional ion-molecule chemistry (see, e.g., 

1995jl. C 2 H 2 is more reliant on efficient ion-molecule 
istry for its formation than HCN because neutral-neutral path¬ 
ways to carbon-chain growth have significant activation barriers 
(see Figure |4). The barriers en route to HCN, in comparison, 
are lower (as discussed in Sect. |3.2| i. Thus, switching off X-ray 
chemistry has a larger effect on the magnitude of the peak abun¬ 
dance reached by C 2 H 2 in the atmosphere as both ion-molecule 
and neutral-neutral pathways are inhibited in the cooler disks. 


Herbst 

chem- 


4.3. On the influence of initial nitrogen reservoirs 


The results have so far suggested that X-ray-induced chemistry is 
crucial for efficient molecular synthesis in the disk atmosphere in 
the planet-forming regions of protoplanetary disks around cool 
stars. An additional scenario to consider is the effect of the as¬ 
sumed initial abundances at the beginning of the calculation. In 
the results presented thus far, a set of initial abundances from 
the output of a dark cloud model were used (see Table |5J. The 
calculations begin with a ratio of N:N 2 :N 2 i ce :NH 3 i ce equal to 
1.0:0.21:0.39:0.26. Taking inspiration from the recent work by 
Schwarz & Bergin (2014) 1, we run an additional set of calcula¬ 
tions in which we assume all species are in atomic form and (i) 
nitrogen is also in atomic form; (ii) nitrogen begins as N 0 gas; 
(iii) nitrogen begins as N 2 ice; and (iv) nitrogen begins as \'H , 
ice. In this way, we investigate the degree of chemical processing 
in the disk atmosphere for each initial nitrogen reservoir. 

In Figure fl4] the fractional abundances of gas-phase N 2 (top 
row), NH 3 (middle row), and HCN (bottom row) are shown as a 
function of Z/R at 1 AU. An equivalent plot for R - 10 AU is pre¬ 
sented in the Appendix (Figure pO) . The results show that the 
calculated abundances in the disk atmosphere (Z/R > 0.1) are in¬ 
dependent of the form of the initial nitrogen reservoir at a radius 
of both 1 and 10 AU for all models. The chemistry in the disk 
atmosphere has achieved steady state by 10 6 years and has ‘for¬ 
gotten’ its origins. On the other hand, the abundances in the disk 
midplane are very sensitive to the initial nitrogen reservoir. The 
results for cases (ii) and (iii) show that similar abundances are 
achieved regardless of whether N, begins in gas or ice form. At 
1 AU the dominant factor is whether nitrogen begins as NH 3 ice 
(case (iv)). Here, the abundance of N 2 in the midplane is lower 
relative to the fiducial model with the difference increasing with 
increasing spectral type. Correspondingly, the abundance of gas- 
phase NH 3 is significantly higher increasing by around one order 
of magnitude compared with the fiducial model. At the other ex¬ 
treme, beginning with all species in atomic form generates the 
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lowest abundance of NH 3 in the disk midplane, especially in the 
case of the M dwarf and T Tauri disks with differences between 
1 and 6 orders of magnitude when compared with the model in 
which nitrogen begins as ammonia ice. NH 3 ice is thought to be 
produced in situ on the surfaces of grains or on or within the ice 
mantle; hence, the efficiency of the conversion from atomic N to 
NH 3 is very sensitive to temperature and a sufficient flux of both 
atomic N and H must reside on the surface for the reaction to 
proceed. 

The initial nitrogen reservoir also has an effect on the mid¬ 
plane HCN abundance. For the T Tauri and Herbig Ae disks, at 
1 AU, HCN reaches the highest abundance for the atomic model, 
followed by the model in which N begins as ammonia ice. Both 
models which have nitrogen initially in the form of N 2 produce 
the lowest abundance of HCN in the disk midplane. The results 
at 10 AU show less spread except in the case of the Herbig Ae 
disk. NH 3 and HCN are both significantly enhanced in the mid¬ 
plane for the model in which N begins as NH 3 ice. 

At elevated temperatures, HCN is formed via the reaction 
between H-, and CN (|Baulch et al.|1994|), the latter of which has 
various routes to formation viaatomic nitrogen and nitrogen hy¬ 
drides, e.g., N + CH or C + NH. Both examples are barrierless 
reactions (|Brownsword et al.| 1 996| [Smith et al.|2004[ |Daranlot| 


et al.|2013) . Conversely, reactions which directly produce CN or 


HCN fromN 2 (e.g., C + N 2 or CH + N 2 ) have large reaction bar- 


riers (» 10,000 K, Baulch et al .||1994[ Rodgers & Smith|1996| >. 
In the dense, warm midplane, NH 3 is more easily broken apart 
by cosmic-ray-induced photodissociati on than N 2 thus re leasing 
nitrogen hydrides into the gas phase (Gred el et al.|1989||Heays| 
|et al.|2014 1 . There is also a direct (and barrierless) route to HCN 
from NH,, 


NH 3 + CN 


HCN + NH-, 


(4) 


(Sims et al. 1994 1 . Hence, the nitrogen chemistry in the disk 


midplane (and resultant abundance and distribution of N-bearing 
molecules) is sensitive to whether nitrogen begins in the form of 
N 2 or NH 3 with the latter resulting in an increase in N-bearing 
species in the midplane. 


4.4. Comparison with previous models 


Astrochemical models of the inner regions of protoplanetary 
disks have often been neglected in favour of the outer disk 
(>> 10 AU), motivated by the larger molecular inventory ob¬ 
served via emission at (sub)mm wavelengths (see, e.g., the re¬ 
cent review by Dutrey et af|[ 2014j i. However, focus will return 
to the planet-forming regions of protoplanetary disks driven by 
spatially-resolved observations of molecules at (sub)mm wave¬ 
lengths with the ALMA Large Millimeter/Submillimeter Array 
(ALMA, see, e.g., Qi et al. 2013[) and the launch of the James 


Webb Space Telescope (JWST) in 2018 (see, e.g., Gardner et al. 
[2006) . 


Early models of the planet-forming region spanned a wide 
range of complexity in both physics and chemistry and have pri¬ 
marily focussed on disks around T Tauri stars. [Willacy et al.[ 
( 1998j > used a one-dimensional dynamical model to determine 
the chemical composition of the midplane from 0.1 to 100 AU. 
They used a chemical network which included gas-phase chem¬ 
istry and gas-grain interactions (freezeout and desorption). They 
conclude that neutral-neutral chemistry is more important than 
ion-molecule chemistry for controlling the abundances in the 
midplane, similar to that found here. They also concluded that 
the chemical composition of the inner midplane was not depen¬ 


dent on the initial molecular abundances. This is a different con¬ 
clusion to that presented here where we find that the composition 
of the initial ice reservoir plays a crucial role in the subsequent 
gas-phase chemistry; however, the model presented here uses 
a chemical network which is significantly more expansive than 
that adopted in Willacy et al.| ( |1998] l and includes grain-surface 
chemistry and ice mantle processing. 

|Markwick et al.| ( |2Q02] ) calculated the two-dimensional 
chemical composition of the inner region assuming a disk heated 
internally by viscous heating only, and using a similar chemical 
network to that from |Willacy et al 7| ( ] 1998) . They also included a 
simple prescription for the X-ray ionisation rate throughout the 
vertical extent. In general, |Markwick et al. ( |2Q02[ > compute sig¬ 
nificantly higher column densities for C 2 H 2 than found in this 
work and they also do not produce OH in the disk surface layer. 
This is because they neglected photodissociation and photoioni¬ 
sation by photons originating from the central star and the exter¬ 
nal interstellar radiation field, processes that we find are impor¬ 
tant for governing the chemistry in the ‘observable’ molecular 
layer in the inner regions. 

More recently, Agundez et al.| ( |2008 i explored the chem¬ 
istry in the inner regions of a disk around a T Tauri star using 
a model similar to that adopted for photon-dominated regions 
(PDRs). They find similar conclusions to here: temperature- 
activated neutral-neutral chemistry helps to build chemical com¬ 
plexity in the disk atmosphere. However, we also find that X-ray 
driven chemistry is potentially very important for building addi¬ 
tional complexity, by releasing free atomic (or ionic) carbon and 
nitrogen into the gas deeper into the disk atmosphere. Qualita¬ 
tively, we see the same behaviour in the column densities in the 
disk atmosphere with radius: the column density starts higher 
then decreases at a radius specific to each molecule. C,H 2 de¬ 
creases first followed by HCN, then H 0 O. OH, on the other hand 
remains flat. The peak column densities calculated for C 2 H 2 and 
HCN by Agundez et al. (2008 1 (~ 10 16 cm -2 ) are somewhat 
lower than those computed here for the T Tauri disk. We also 
find much larger column densities for both OH and H 2 0 in the 
atmosphere. The reason for the particularly low OH column den¬ 
sity is unclear but Agundez et al.| ( |2008) > do neglect heating via 
UV excess emission from the star and also assume that the gas 
and dust temperatures are equal; whereas we find that the gas is 
significantly hotter than the dust in the region of the inner disk at¬ 
mosphere where the molecules reside. The column densities that 
we calculate for OH and H 2 0 for the T Tauri disk are also in line 
with those determined in the work by Bethell & Bergin (2009) 
and Glassgold et al. (|2009[); however, unlike Bethell & Bergin 


(2009), we do not consider dust grain settling and so do not need 
to invoke H 0 O self shielding as a mechanism for explaining the 
survival of gas-phase water in the disk atmosphere. It is interest¬ 
ing that Bethell & Bergin| ( |2009| l (and follow up work by |Du &| 
Bergin|2014 1 also predict an increase in column density of H 2 0 


in the disk atmosphere with increasing stellar FUV luminosity 
(and corresponding hardening of the radiation field). 


4.5. A connection between the disk atmosphere and the 
planet-forming midplane? 

Currently, the only means to probe the inner planet-forming re¬ 
gions of protoplanetary disks is via near- to mid-infrared obser¬ 
vations. Future high-angular-resolution observations at (sub)mm 
wavelengths with ALMA may also elucidate the molecular com¬ 
position of the inner region; however, due to the high column 
densities in the inner disk, dust opacity begins to affect even 

Article number, page 15 of|28| 











































































A&A proofs: manuscript no. inner_disk_chemistry_accepted 


M Dwarf Disk T Tauri Disk 







- No 

- No 

- - ■ Fid 

—-- - 

shielding hCN 

X-rays 

ucial 

y 

/X\ 


Herbig Ae Disk 




0.0 0.1 0.2 0.3 0.4 


Z/R 


Fig. 12. Fractional abundance with respect to gas number density of N (top row), N 2 (middle row), and HCN (bottom row) as a function of Z/R at 
R = 1 AU for each disk model. The black dashed lines, purple solid lines, and green solid lines represent results from the fiducial model (including 
N-, shielding and X-rays), the model with N-> switched off, and the model with X-rays switched off, respectively. 
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Fig. 13. Fractional abundance with respect to gas number density of C 2 FF as a function of Z/R at R = 1 AU for the M dwarf and T Tauri disk 
models. The black dashed lines and orange solid lines represent results from the fiducial model (including N-, shielding and X-rays) and the model 
with X-rays switched off, respectively. 


(sub)mm line emission. The higher sensitivity and spectral res¬ 
olution of JWST will also allow the measurement of absorption 
features by gases and ices other than water in nearby edge-on 
protoplanetary disk systems. 


Given that forming planets sweep up material primarily from 
the disk midplane, it is worth exploring to what degree the 
composition of the disk atmosphere reflects that of the mid¬ 
plane. During the planet-formation process, the molecules ac¬ 
creted (whether gas or ice) are ultimately reprocessed in the 
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Fig. 14. Fractional abundance with respect to gas number density of N 2 (top row), NH 3 (middle row), and HCN (bottom row) as a function of Z/R 
at R = 1 AU for each disk model. The black dashed lines, orange lines, cyan lines, green lines, and gold lines represent results from the fiducial 
model and a model in which (i) all species are initially atomic; (ii) all nitrogen is in N 2 gas; (iii) all nitrogen is in N, ice; and (iv) all nitrogen is in 
NH 3 ice, respectively. 


planet atmosphere. Recent population synthesis models suggest 
that the main contribution to the heavy element content in the 
atmospheres of forming planets are ices accreted during the for¬ 
mation of the planetary embryo and icy planetesimals which are 
captured during the gas accretion stage and vaporised in the at¬ 
mosphere (see, e.g., |Thiabaud et al.|2014] >. 

In Figure[l5] the percentage contribution of each of the dom¬ 
inant molecular carriers of oxygen (top row), carbon (middle 
row), and nitrogen (bottom row) are plotted as a function of ra¬ 
dius. In Figure [16] we present the equivalent data for the disk 
atmosphere only. The dominant carriers are identified as those 
species which contribute most to the total column density of each 
element (for the data plotted in Figure[T5j) and to the column den¬ 
sity in the disk atmosphere, down to the r(14 fim) = 1 surface 
(for the case of the data plotted in Figure [l6]). The data identified 
as ‘Other’ refers to the summation over all other species which 
contribute to the elemental abundance, but which individually 
do not contribute significantly. The identity of ‘Other’ depends 
on the location in the disk and the particular disk model but are 
typically complex organic ices (e.g., H 2 CO, CH 3 OH, etc.) in the 
outer regions of the two cooler disks, and gas-phase hydrocar¬ 
bons (e.g., C n H m ) in the inner regions of the two warmer disks. 
For the most part, the plotted species contribute * 90 - 100% to 
the total elemental abundance. 


The dominant carriers throughout the vertical extent is set 
by the locations of snow lines: gas-phase H,0 in the inner re¬ 
gions is superseded by H 2 0 ice once the snow line is surpassed. 
An interesting result for the M dwarf disk is that CO, ice car¬ 
ries most of the oxygen beyond ~ 3 AU indicating an efficient 
conversion from H,0 ice to CO, ice on and within ice man¬ 
tles. This is not seen in the two warmer disks because the higher 
temperatures make it difficult for CO to reside on the grains suf¬ 
ficiently long for reaction with OH. Beyond the H,0 snowline, 
the dominant gas-phase carrier of elemental oxygen is CO: thus 
the gas-phase C/O ratio seemingly tends towards 1, whereas the 
ice mantle remains oxygen-rich. However, in each disk, there 
are regions where molecular oxygen has a non-negligible con¬ 
tribution to the total oxygen column, on the order of 10%. The 
two-dimensional fractional abundances of important gas-phase 
oxygen-bearing molecules not discussed in detail in the text. 


such as CO, 0 2 , and CO,, are presented in Figure A.3 in the 
Appendix. 

The dominant carbon carrier is CO; in the outer disk of the 
cooler M dwarf, CO, ice and other complex organic ices such 
as CH 3 OH take over. In the inner region of the T Tauri and Her¬ 
big Ae disks, gas-phase CH 4 and other hydrocarbons begin to 
contribute at the level of 30 - 40%: the physical conditions in 
the inner regions of the hotter disks are such that the chemistry 
tends towards thermochemical equilibrium in which hydrocar- 
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bons dominate over CO. In the outer regions of the two hotter 
disks, again either C0 2 ice or gas also contributes at the level of 


30%. For the disk around TW Hya, Bergin et al. ( [2014 ) suggest 
that gas-phase CO is an order of magnitude lower in abundance 
than that expected if 100% of the freely available (i.e., not con¬ 
tained in refractory material) elemental carbon were locked up 
in gas-phase CO. Our results for the M dwarf disk show a de¬ 
pletion in CO in the outermost regions consistent with this level 
of depletion; however, we find that significantly more carbon is 
contained in CO gas in the warmer disks (~ 40 - 50%, at least 
within a radius of 10 AU). 

The picture for nitrogen is more simple: in all cases, gas- 
phase N ? is the primary carrier. In the M dwarf disk, NH 3 con¬ 
tributes at the level of 10% whereas in the two hotter disks, 
there is a region where the contribution from gas-phase NH 3 
approaches 50% which is temperature dependent. This is again 
because the conditions in the inner disk midplane approach 
thermochemical equilibrium. The low abundance of XH, ice 
in this region suggests that planetary atmospheres which gain 
the bulk of their heavy elements from planetesimal accretion 
will be depleted in nitrogen relative to carbon and oxygen. The 
two-dimensional fractional abundances of important gas-phase 
nitrogen-bearing molecules not discussed in detail in the text. 


such as N, and NH 3 , are presented in Figure A.6 in the Ap¬ 


pendix. 

In Figure [16] the equivalent values for the disk atmosphere 
are plotted. For the oxygen carriers, the main difference is the 
increased contribution of gas-phase O, to the oxygen budget in 
the atmosphere (up to « 50%). Where O, dominates the oxygen 
budget depends on the disk model, moving from the innermost 
regions of the M dwarf disk to the outermost region of the Her- 
big Ae disk. O, is formed in the atmosphere via reaction between 
atomic oxygen and OH and destroyed via reaction with atomic 
hydrogen and carbon, with photodissociation increasing in im¬ 
portance as the spectral type of the star increases. Figure[9]shows 
the dominant formation and destruction mechanisms for O, in 
the disk atmosphere, and the two-dimensional abundance dis¬ 
tributions and column densities of atomic and molecular oxygen 
are given in the Appendix. Where O, dominates, CO, also makes 
a non-negligible contribution in the atmosphere (up to ~ 30%). 

The story for carbon is more simple: CO dominates in the 
disk atmosphere over much of the radial extent of all disks, with 
some contribution from CO,, as found for the full column val¬ 
ues. In the two warmer disks, the gas-phase carbon carriers in 
the disk atmosphere switch to ‘Other’, in this case, various gas- 
phase hydrocarbons, C n H m , which is not representative of the 
total column. Gas-phase HCN also has a non-negligible contri¬ 
bution to both the carbon and nitrogen budget in the atmosphere 
(up to 30% and 60% respectively). 

Similarly, gas-phase N 2 dominates the nitrogen budget 
throughout most of the disk atmosphere; however, in the inner¬ 
most regions of the two warmer disks, gas-phase NH 3 and HCN 
take over. Again, this is not representative of the nitrogen bud¬ 
get throughout the disk vertical extent which is dominated by 
N 2 (~ 90%). In Figure[5]we also show those reactions responsi¬ 
ble for the formation ana destruction of N 2 and NH 3 in the disk 
atmosphere. 

In Figure [F7] the C/O ratio is plotted as a function of radius 
for three cases. Also plotted is the assumed underlying elemental 
ratio (C/O = 0.44, gray dashed lines). The left-hand panel shows 
the ratio when summed over the dominant ice reservoirs (H,0 
and C0 2 ice). The ice is more oxygen rich than the underlying 
elemental ratio; however, for the two cooler disks, as CO, begins 
to freezeout and/or form, the C/O ratio of the ice tends towards 


the underlying ratio. On the other hand, the gas is either repre¬ 
sentative of the elemental ratio or more carbon rich (as shown in 
the middle panel). The behaviour of the C/O ratio depends on the 
locations of the H,0 and CO, snow lines which move outwards 
with increasing stellar spectral type. When H,0 ice freezes out 
and/or forms, the C/O ratio in the gas increases as oxygen is re¬ 
moved. Once CO, ice begins to freeze out and/or form, the gas 
C/O ratio begins to decrease again and the ice correspondingly 
becomes relatively more carbon rich. This is most clearly seen 
for the M dwarf disk. In the final panel, the C/O ratio calculated 
using ‘observable’ tracers only is shown (CO, CO,, H,0, C,H 2 , 
CH 4 , and HCN gas). Comparing the final panel with the middle 
panel, for the two cooler disks, the gas in the disk atmosphere ap¬ 
pears significantly more carbon rich when calculated using just 
the listed tracers. Notably, the C/O ratio for the M dwarf disk 
appears > 1 throughout most of the disk atmosphere (> 1 AU). 
This is primarily due to the presence of gas-phase O, in the at¬ 
mosphere which is another ‘hidden’ reservoir of atomic oxygen 
(in addition to H,0 ice and CO, ice). 

The chemical model results suggest that disks around cooler 
stars might appear more carbon rich without the need for addi¬ 
tional sinks (neither chemical nor physical) to account for the 
depletion of oxygen. However, we stress that detailed radia¬ 
tive transfer calculations are necessary to confirm definitively 
whether the chemical models replicate the trends seen in the ob¬ 
servations. Because this is a non-trivial calculation for a two- 
dimensional disk structure requiring careful consideration of the 
dust structure and size distribution, this is beyond the scope of 
the work presented here but is planned future work. Assuming 
that the trends in the chemical models do translate into trends 
in the simulated emission, the derivation of the C/O ratio using 
CO, C0 2 , H,0, CH 4 , C,H,, and HCN alone, may overestimate 
the underlying C/O ratio by up to a factor of two. Similarly, the 
C/N ratio may be overestimated by between one and several or¬ 
ders of magnitude if observations of HCN and NH 3 in the disk 
atmosphere alone are used to determine the underlying C/N ra¬ 
tio. In the case that the dominant source of heavy elements in 
a planetary atmosphere is icy planetesimals rather than the gas, 
then this overestimation increases to a factor of 10. Although less 
extreme than the estimation of the C/N ratio, a factor of two is 
sufficient to incorrectly assume the formation of carbon-rich ver¬ 
sus oxygen-rich planetary atmospheres, and the dominant carbon 
and oxygen carriers subsequently found therein (see, e.g., Mad- 
husudha n et al.|2013| l. 


4.6. On the abundance of 0 2 in protopianetary disks 

The models predict that gas-phase O, is a significant oxygen 
reservoir in the inner regions of protopianetary disks. O, may 
contain up to 10% of the total available oxygen over the full 
vertical extent of the disk, and this increases to up to 50% when 
the disk atmosphere only is considered (see Figures [T5| and [T6]>. 

O, has proven to be a somewhat elusive molecule in inter¬ 
stellar and circumstellar environments. Upper limits determined 
by SWAS (Submillimeter Wave Astronomy Satellite) and Odin 


towards nearby molecular clouds constrain 0 2 /H 2 < 10 '(Gold¬ 


smith et ak| |2000[ Pagani et al.| |2003| >. On the other hand, ISO 
upper limits for the abundance of 0 2 ice in dark clouds are 
much more conservative (0,/H,0< 0.6, Vandenbussche et al. 
1999jl. However, gas-phase 0 2 has been successfully detected 
towards two warmer sources, the dense core, p Oph A (Lars- 
son et al.|2007[ |Liseau et al.|2012]>, and Orion ( Goldsmith et al. 
201 1| >. Furthermore, Yildiz et al 020131 report a deep Herschel- 
determined upper limit towards the low-mass Class 0 protostar. 
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Fig. 15. Percentage contribution of the dominant molecular carriers to the total elemental oxygen, carbon, and nitrogen column densities, as a 
function of radius. 


NGC 1333-IRAS 4A (0 2 /H 2 < 6 x 10 9 ). The authors find that 
Ot is absent in both the outer cold envelope and inner hot core 
and conclude that the material entering protoplanetary disks is 
likely poor in molecular oxygen (gas and ice). 

Early gas-phase only chemical models routinely overpre¬ 
dicted the abundance of 0 2 in dark clouds, which was postulated 
to form primarily via ion-molecule chemistry (see, e.g., |Bergtn| 
et aT| 2000) . More modern and sophisticated gas-grain models 
are able to reproduce the low gas-phase abundance of Oi pro¬ 
vided the conversion of 0 2 ice into H 2 0 ice is included, and the 
chemistry is allowed to evolve for sufficiently long time scales 
(see, e.g., |Bergin et al. 12000] |Roberts & Herbst|2002[ Yildiz et 


|al.|2013| ). Recent laboratory experiments have shown that the hy 
drogenation of 0 2 ice is rapid at low temperatures (loppolo et al. 


2008} Miyauchi et al.[2008 l. 


The origin of gas-phase 0 2 in the inner regions of protoplan¬ 
etary disks is different to that expected in dark clouds. In both 
cases, the vital reaction is formation via O + OH. This reaction 
has been well studiecQ across the temperature range of interest 
for circumstellar environments. However, the origin of OH in 
dark clouds is via the dissociative recombination of H 3 0 + which 
is generated via successive proton-donation reactions originat¬ 
ing from O + H 3 . In the inner regions of disks, OH is gener¬ 
ated via the reaction between O + H 2 , and atomic oxygen, in 
turn is released from CO via photodissociation. X-ray-induced 
dissociation, or via reaction with He + . 0 2 is destroyed via pho¬ 
todissociation and reactions with C and H to yield CO and OH, 


respectively (see Figure[9]). These latter two reactions have well- 
constrained rate coefficients (Gep pert et al.||2000) |Miller et al.| 
2005). Gas-phase 0 2 is able to persist in the disk atmosphere for 


the same reasons as gas-phase HiO: the gas temperature is suf¬ 
ficiently high to activate the required neutral-neutral chemical 
reactions. The model T Tauri disk generates column densities of 
OH and H 2 0 in good agreement with those observed towards T 
Tauri stars (see Figure [TO}: hence, given that the rate coefficients 
for the important reactions are well studied, the prediction that 
gas-phase O-, may also be relatively abundant in T Tauri disks is 
substantive. 


5. Summary 

In this work, the chemistry and resulting molecular composi¬ 
tion of the planet-forming regions (< 10 AU) of protoplanetary 
disks has been explored, with the aim to investigate potential 
reasons for the trends seen in near- to mid-infrared observations. 
The results demonstrate that, as the effective temperature of the 
central star increases, the molecular complexity of the disk at¬ 
mosphere decreases, showing that the FUV luminosity of the 
host star plays an important role in determining the molecular 
composition of the disk atmosphere. The weaker FUV flux im¬ 
pinging upon disks hosted by M dwarf stars allows molecules 
to reach relatively high abundances in the atmosphere: X-ray- 

5 http://kida.obs.u-bordeauxl.fr/ 
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Fig. 16. Same as Figure[l5]for the disk atmosphere only. 
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Fig. 17. C/O ratio as a function of radius for the total ice column density (left panel), the total gas column density (middle panel), and for the 
‘observable’ disk atmosphere only (right panel). The gray dashed line indicates the underlying C/O elemental ratio (0.44). 


induced chemistry can further increase molecular complexity by 
driving a rich ion-molecule chemistry, helping to qualitatively 
explain the large column densities of small organic molecules, 
such as, CtH 0 and HCN seen in M dwarf disk atmospheres. The 
key process is the liberation of free carbon and nitrogen from 
their main molecular reservoirs (CO and N,„ respectively) via 
X-rays. The results shown here, in conjunction with the mid- 
IR observations, suggest that these objects are good targets for 
ALMA which can further help elucidate chemistry in disks to¬ 
wards the low-mass low-luminosity regime. 


The chemical models suggest that the gas in the inner re¬ 
gions of M dwarf disks is generally more carbon rich than that 
in disks around T Tauri stars. When the C/O gas-phase ratio is 
calculated using only observable tracers in the disk atmosphere, 
then the ratio appears larger than it actually is (and C/O —> 1). 
This is because the models predict that gas-phase 0 2 is a sig¬ 
nificant reservoir of oxygen in the disk atmosphere beyond the 
water snowline. 

The results also demonstrate a degree of chemical decou¬ 
pling between the disk atmosphere and the midplane. The gas 
is generally more carbon-rich than the midplane ices (see Fig- 
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ure E5- This is further corroborated by our studies on the impor¬ 
tance of the initial nitrogen reservoir (whether atomic nitrogen, 
molecular nitrogen, or ammonia). We find that this does not play 
a role in determining the resulting composition of the observable 
molecular layer: the chemistry is at steady state. However, the 
initial reservoir is important for determining the composition in 
the disk midplane where, generally, the chemical timescales are 
longer. For example, icy planetesimals forming in disks where 
much of the initial nitrogen is locked up in \'H ( ice may be 
more nitrogen rich than those for which nitrogen was contained 


primarily in the more volatile N, (see also Schwarz & Bergin 
|20l4) . 

Whether the trends seen in the chemical models can repro¬ 
duce those derived from observations remains to be confirmed 
via calculations of the molecular emission and will be conducted 
in future work. Assuming that the chemical trends do translate 
into observable trends, near- to mid-IR observations of the dom¬ 
inant tracers in the atmosphere may overestimate the underlying 
C/O and C/N ratios of the gas and ice in the region in which 
forming planets sweep up most of their material by up to a fac¬ 
tor of 10 and more than an order of magnitude, respectively. 

Despite the results qualitatively demonstrating some of the 
observed trends, several issues remain. The large column densi¬ 
ties of water vapour and correspondingly strong water emission 
lines at near- to mid-IR wavelengths predicted by models of Her- 
big Ae disks have not been confirmed by observations, indicating 
that there may be a heretofore unconsidered destruction mech¬ 
anism for gas-phase water at high temperatures. For example, 
reactions of rovibrationally excited OH with atomic hydrogen 
may shift the ratio of O/OH/H.O in the atmosphere. Alternative 
hypotheses are that the molecular line emission is veiled by the 
strong stellar continuum emission (as discussed in Pontoppidan 


et al.|2010|> or that the disk dust structure plays an important role. 


A further issue is that the M dwarf disk model predicts a 
C 2 H 2 /HCN ratio which is an order of magnitude lower relative 
to the observations. An increase in the underlying C/O or C/N 
elemental ratio in the disk atmosphere may help explain the high 
C 2 H 2 /HCN ratio in M dwarf disks with the enrichment in carbon 
relative to oxygen and nitrogen caused by vertical or radial mix¬ 
ing. In this scenario, less volatile species (e.g., H z O and NH 3 ) are 
transported to regions where it is sufficiently cold for freezeout 
onto dust grains (see, e.g., Stevenson & Luninep~988|; Meijerink 


|et al.||2009] >. If the ice-coated dust grains are sufficiently large, 
they become decoupled from the gas and settle to the disk mid¬ 
plane. In this way, the gas can become enriched in more volatile 
species, e.g., CO, which alters the underlying elemental balance 
of the atmosphere. However, as shown in the results here, molec¬ 
ular nitrogen (which is volatile) is significantly more abundant 
in the disk atmosphere than NH 3 . Hence, this mechanism is un¬ 
likely to enrich the disk atmosphere in carbon relative to nitro¬ 
gen but perhaps the chemistry of nitrogen-bearing species is per¬ 
turbed by the depletion of oxygen via this mechanism. 

The gas-phase chemical network used for HCN and C 2 H 2 
is (to our knowledge) relatively complete, with experimentally 
measured and/or critically reviewed reaction rate coefficients 
used where available. We have shown that self-shielding of N 2 
alone is not sufficient to explain C 2 H 7 /HCN in the inner regions 
and has a more significant effect in disks around stars with a 
higher FUV luminosity. It remains to be confirmed whether this 
also holds for disks in which grain growth and settling have gen¬ 
erated a relatively dust-poor atmosphere for which molecular 
shielding dominates over dust shielding. Instead, X-ray-induced 
chemistry is more important for releasing atomic nitrogen from 
Nt in the M dwarf disk. Furthermore, we find that the exclu¬ 


sion of X-ray-induced chemistry only increases the discrepancy 
with observation, because the formation of C 2 H 2 is primarily 
driven by ion-molecule chemistry in the M dwarf disk atmo¬ 
sphere, as shown in Figures [4] and |T3] However, all chemical 
networks suffer from a degree of uncertainty; hence, a systematic 
sensitivity study of the relative abundance of these two species 
over the parameter space of physical conditions in protoplan¬ 
etary disks around cool stars, is worthy of further exploration. 
This would also confirm the hypothesis presented here, that is, 
that X-ray chemistry is responsible for the higher C-,H 2 /HCN 
ratio observed in cool stars. 

We finish by stating that the future is bright for near- to mid- 
IR astronomy with MIRI (Mid-InfraRed Instrument, Wright et 
al.|2004[ ) on IWST and METIS ( [Brandi et al.|20l4| , currently be¬ 


ing developed for installation on the European Extremely Large 
Telescope (E-ELT). Both facilities will have significantly higher 
spectral resolution (~ 3000 and 100,000 respectively) across the 
wavelength range of interest for probing the molecular compo¬ 
sition of the planet-forming regions of nearby protoplanetary 
disks. 
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Appendix A: Supplementary figures 

In Figures |AJ~| and [A~2j we show results at a radius of 10 AU 
from the same models described in Figures [12] and [T4| 

In Figures [A3] to |A.6~1 we show the fractional abundances as 
a function of radius, R, and height scaled by the radius, Z/R , for 
those species not discussed in detail in the text. These include 
important oxygen- (e.g., 0 2 , CO, and CO-,), carbon- (e.g., C n H m ) 
and nitrogen-bearing species (e.g., N 2 and NH 3 ). 



n(X)/n n(X)/n n(X)/n 


Walsh, C. et al.: The molecular composition of protoplanetary disks 


M Dwarf Disk 


T Tauri Disk 


Herbig Ae Disk 


1 


1 


1 


10' 3 

10' 4 

10' 5 

10' 6 

10" 7 

10' 8 

10' 9 

O ' 10 

10' 3 

10' 4 

10' 5 

10' 6 

10' 7 

10' 8 

10' 9 

O ' 10 

10' 3 

10' 4 

10' 5 

10' 6 

10- 7 

10' 8 

10' 9 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


Z/R 


Z/R 


Z/R 


Fig. A.l. Fractional abundance with respect to gas number density of N (top row). N 0 (middle row), and F1CN (bottom row) as a function of 
Z/R at R = 10 AU for each disk model. The black dashed lines, purple solid lines, and green solid lines represent results from the fiducial model 
(including N, shielding and X-rays), the model with N, switched off, and the model with X-rays switched off, respectively. 
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Fig. A.2. Fractional abundance with respect to gas number density of N-, (top row), NH 3 (middle row), and HCN (bottom row) as a function of Z/R 
at R = 10 AU for each disk model. The black dashed lines, orange lines, cyan lines, green lines, and gold lines represent results from the fiducial 
model and a model in which (i) all species are initially atomic; (ii) all nitrogen is in N, gas; (iii) all nitrogen is in N, ice; and (iv) all nitrogen is in 
NH 3 ice, respectively. 


Article number, page 24 of 


28 

























Walsh, C. et al.: The molecular composition of protoplanetary disks 


M Dwarf Disk 


T Tauri Disk 


Herbig Ae Disk 




: so 2 







R (AU) R (AU) 



Fig. A.3. Fractional abundance relative to total gas number density of oxygen-bearing species for the M dwarf disk (left-hand column), T Tauri 
disk (middle column), and Herbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density (integrated from 
the surface downwards) at which r » 1 at 3 fim and 14 jum, respectively. 
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Fig. A.4. Fractional abundance relative to total gas number density of carbon-bearing species for the M dwarf disk (left-hand column), T Tauri 
disk (middle column), and Herbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density (integrated from 
the surface downwards) at which r » 1 at 3 pm and 14 pm, respectively. 
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Fig. A.5. Fractional abundance relative to total gas number density of carbon-bearing species for the M dwarf disk (left-hand column), T Tauri 
disk (middle column), and Herbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density (integrated from 
the surface downwards) at which r » 1 at 3 fim and 14 /rm, respectively. 
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Fig. A.6. Fractional abundance relative to total gas number density of nitrogen-bearing species for the M dwarf disk (left-hand column), T Tauri 
disk (middle column), and Flerbig Ae disk (right-hand column). The dotted and dot-dashed lines indicate the dust column density (integrated from 
the surface downwards) at which r » 1 at 3 pm and 14 pm, respectively. 
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